SIX- VERTEX, LOOP AND TILING MODELS: 
INTEGRABILITY AND COMBINATORICS 
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Abstract. This is a review (including some background material) of the author's work and related 
activity on certain exactly solvable statistical models in two dimensions, including the six-vertex 
model, loop models and lozenge tilings. Applications to enumerative combinatorics and to algebraic 
geometry are described. 
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Introduction 

Exactly solvable (integrable) two-dimensional lattice statistical models have played an important 
role in theoretical physics: starting with Onsager's solution of the Ising model, they have provided 
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non-trivial examples of critical phenomena in two dimensions and given examples of lattice real- 
izations of various known conformal field theories and of their perturbations. In all these physical 
applications, one is interested in the thermodynamic limit where the size of the system tends to 
infinity and where details of the lattice become irrelevant. In the most basic scenario, one considers 
the infra-red limit and recovers this way conformal invariance. 

On the other hand, combinatorics is the study of discrete structures in mathematics. Combina- 
torial properties on integrable models will thus be uncovered by taking a different point of view 
on them which considers them on finite lattices and emphasizes their discrete properties. The 
purpose of this text is to show that the same methods and concepts of quantum integrability lead 
to non-trivial combinatorial results. The latter may be of intrinsic mathematical interest, in some 
cases proving, reproving, extending statements found in the literature. They may also lead back to 
physics by taking appropriate scaling limits. 

Let us be more specific on the kind of applications we have in mind. First and foremost comes the 
connection to enumerative combinatorics. For us the story begins in 1996, when Kuperberg showed 
how to enumerate alternating sign matrices using Izergin's formula for the six-vertex model. The 
observation that alternating sign matrices are nothing but configurations of the six- vertex model in 
disguise, paved the way to a fruitful interaction between two subjects which were disjoint until then: 
(i) the study of alternating sign matrices, which began in the early eighties after their definition by 
Mills, Robbins, and Rumsey in relation to Dodgson condensation, and whose enumerative properties 
were studied in the following years, displaying remarkable connections with a much older class of 
combinatorial objects, namely plane partitions; and (ii) the study of the six-vertex model, one of 
the most fundamental solvable statistical models in two-dimensions, which was undertaken in the 
sixties and has remained at the center of the activity around quantum integrable models ever since. 

One of the most noteworthy recent chapters in this continuing story is the Razumov-Stroganov 
conjecture, in 2001, which emerged out of a collective effort by combinatorialists and physicists 
to understand the connection between the aforementioned objects and another class of statistical 
models, namely loop models. The work of the author was mostly a byproduct of various attempts 
to understand (and possibly prove) this conjecture. A large part of this manuscript is dedicated to 
reviewing these questions. 

Another interesting, related application is to algebraic combinatorics, due to the appearance of 
certain families of polynomials in quantum integrable models. In the context of the Razumov- 
Stroganov, they were introduced by Di Francesco and Zinn-Justin in 2004, but their true meaning 
was only clarified subsequently by Pasquier, creating a connection to representation theory of affine 
Hecke algebras and to previously studied classes of polynomials such as Macdonald polynomials. 
These polynomials satisfy relations which are typically studied in algebraic combinatorics, e.g. 
involving divided difference operators. The use of specific bases of spaces of polynomials, which 
is necessary for a combinatorial interpretation, connects to the theory of canonical bases and the 
work of Kazhdan and Lusztig. 

Finally, an exciting and fairly new aspect in this study of integrable models is to try to find an 
algebro-geometric interpretation of some of the objects and of the relations that satisfy. The most 
naive version of it would be to relate the integer numbers that appear in our models to problems 
in enumerative geometry, so that they become intersection numbers for certain algebraic varieties. 
A more sophisticated version involves equivariant cohomology or K-theory, which typically leads to 
polynomials instead of integers. The connection between integrable models and certain classes of 
polynomials with geometric meaning is not entirely new, and the work that will be described here 
bears some resemblance, as will be reminded here, to that of Fomin and Kirillov on Schubert and 
Grothendieck polynomials. However there are also novelties, including the use of the multidegree 
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technology of Knutson et al, and we apply these ideas to a broad class of models, resulting in new 
formulas for known algebraic varieties such as orbital varieties and the commuting variety, as well 
as in the discovery of new geometric objects, such as the Brauer loop scheme. 

The presentation that follows, though based on the articles of the author, is meant to be essen- 
tially self-contained. It is aimed at researchers and graduate students in mathematical physics or 
in combinatorics with an interest in exactly solvable statistical models. For simplicity, the inte- 
grable models that are defined are based on the underlying affine quantum group Uq{5l{2)), with 
the notable exception of the discussion of the Brauer loop model in the last section. Furthermore, 
only the spin 1/2 representation and periodic boundary conditions are considered. There are in- 
teresting generalizations to higher rank, higher spin and to other boundary conditions of some of 
these results, on which the author has worked, but for these the reader is referred to the literature. 

The plan of this manuscript is the following. In section [U we discuss free fermionic methods. 
Though free fermions in two dimensions may seem like an excessively simple physical model, they 
already provide a wealth of combinatorial formulae. In fact they have become extremely popular 
in the recent mathematical literature. We shall apply the basic formalism of free fermions to in- 
troduce Schur functions, and then spend some time reviewing the properties of the latter, because 
they will reappear many times in our discussion. We shall then briefly discuss the application to 
the enumeration of plane partitions. Section [2] covers the six-vertex model, and in particular the 
six- vertex model with domain wall boundary conditions. We shall discuss its quantum integrability, 
which is the root of its exact solvability. Then we shall apply it to the enumeration of alternating 
sign matrices. In section [3l we shall discuss statistical models of loops, their interrelations with the 
six-vertex model, their combinatorial properties and formulate the Razumov-Stroganov conjecture. 
Section m introduces the quantum Knizhnik-Zamolodchikov equation, which will be used to recon- 
nect some of the objects discussed previously. The last section, [5l will be devoted to a brief review 
of the current status on the geometric reinterpretation of some of the concepts above, focusing on 
the central role of the quantum Knizhnik-Zamolodchikov equation. 

1. Free fermionic methods 

As mentioned above, we want to spend some time defining a typical free fermionic model and to 
apply it to rederive some useful formulae for Schur functions, which will be needed later. We shall 
also need some formulae concerning the enumeration of plane partitions, which will appear at the 
end of this section. 

1.1. Definitions. 

1.1.1. Operators and Fock space. Consider a fermionic operator ip{z): 

(1.1) v(^)= Yl V'-fc^'"^. ^*i^)= E ^fc^'"' 

with anti-commutation relations 

(1-2) irr, i'sU = Srs [A, i^s]+ = irr, Ash = 

^^{z) and ip*{z) should be thought of as generating series for the ipk a-nd so that z is just 
a formal variable. What we have here is a complex (charged) fermion, with particles, and anti- 
particles which can be identified with holes in the Dirac sea. These fermions are one-dimensional, 
in the sense that their states are indexed by (half-odd-)integers; ip'^ creates a particle (or destroys 
a hole) at location k, whereas ipk destroys a particle (creates a hole) at location k. 
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We shall explicitly build the Fock space T and the representation of the fermionic operators now. 
Start from a vacuum |0) which satisfies 

(1.3) V'fc|0)=0 A;>0, V'fc|0) = A; < 

that is, it is a Dirac sea filled up to location 0: 



|0) 



Then any state can be built by action of the ■0fc and ■0^ from 1 0) . In particular one can define 
more general vacua at level £ G Z: 



(1.4) 



2 2 



•••V'l |o) 
• |o) 



^ > 
^ < 



which will be useful in what follows. They satisfy 

(1.5) Vfc 1^) = k>l, V'fc 1^) = k<i 

More generally, define a •partition to be a weakly decreasing finite sequence of non-negative integers: 
Ai > A2 > • • • > A.„ > 0. We usually represent partitions as Young diagrams (also called Ferrers 
diagrams): for example A = (5, 2, 1, 1) is depicted as 



A 



To each partition A 
(1.6) 



(Ai, . . . , A„) we associate the following state in J^f. 



n) 



Note the important property that if one "pads" a partition with extra zeroes, then the corresponding 
state remains unchanged. In particular for the empty diagram 0, \0;i) = For i = we just 
write |A;0) = |A). 

This definition has the following nice graphical interpretation: the state | A; £) can be described by 
numbering the edges of the boundary of the Young diagram, in such a way that the main diagonal 
passes between i — ^ and £+^', then the occupied (resp. empty) sites correspond to vertical (resp. 
horizontal) edges. With the example above and ^ = 0, we find (only the occupied sites are numbered 
for clarity) 



o-l-o-l-o? 



'-3 
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The where A runs over all possible partitions (two partitions being identified if they are 

obtained from each other by adding or removing zero parts), form an orthonormal basis of a 
subspace of J- which we denote by J-i. ip^ and are Hermitean conjugate of each other. 
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Note that (jl.6p fixes our sign convention of the states. In particular, this implies that when one 
acts with ip). (resp. ip^) on a state |A) with a particle (resp. a hole) at k, one produces a new state 
I A') with the particle removed (resp. added) at k times —1 to the power the number of particles to 
the right of k. 

The states A can also be produced from the vacuum by acting with -0 to create holes; paying 
attention to the sign issue, we find 

(1-7) |A; i) = (-1)I'IV',_a;+i • • • ^i-x'^+m-i K + m) 

where the A^ are the lengths of the columns of A, |A| is the number of boxes of A and m = Ai. 
This formula is formally identical to p.6|) if we renumber the states from right to left, exchange 
ip and and replace A with its transpose diagram A' (this property is graphically clear). So the 
particle-hole duality translates into transposition of Young diagrams. 

Finally, introduce the normal ordering with respect to the vacuum |0): 



(1.8) :^jJ*^p^: =-:^^^^: = 

which allows to get rid of trivial infinite quantities. 



^p*iJk j > 



1.1.2. g[(oo) and £|[(1) action. The operators tp*{z)ip{w) give rise to the Schwinger representation 
of 0l(oo) on whose usual basis is the '■ipr'^s : , t, s G Z + i, and the identity. In the first quantized 
picture this representation is simply the natural action of 0l(oo) on the one-particle Hilbert space 
C^^2 and exterior products thereof. The electric charge Jq = :'ip*'ipr- is a conserved number 
and classifies the irreducible representations of 0l(oo) inside J-', which are all isomorphic. The 
highest weight vectors are precisely our vacua i G Z, so that T = with Ti the subspace 

in which Jq= I. 

The g[(l) current 
(1.9) 3{z)= ■.r{z)^{zy. =Y,Jn 



-n-l 



with Jn = Xlr '■'4'r-n'^r'- forms a 0[(1) (Heiscnberg) sub-algebra of 0[(oo): 

(1.10) [Jm,Jn]=m5n 



-'m,—n 



Note that positive modes commute among themselves. This allows to define the general "Hamil- 
tonian" 

oo 

(1.11) H[t] = ^t,J, 

q=l 

where t = (ti, . . . , tg, . . .) is a set of parameters ("times"). 

The Jq, q > 0, displace one of the fermions q steps to the left. This is expressed by the formulae 
describing the time evolution of the fermionic fields: 

e^WV(^)e-^t*l =e~^'°°=^*^'V(^) 
^^■^^^ e^W^*(z)e-^W = e+^r^i*.-V*(^) 

(proof: compute [Jq,%p^*\z)] = ibz'^^W(z) and exponentiate). Of course, similarly, J_q, q > 0, 
moves one fermion q steps to the right. 
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1.2. Schur functions. 

1.2.1. Free fermionic definition. It is known that the map 1$) ^ (^|e^W 1$) 

is an isomorphism 

from J^i to the space of polynomials in an infinite number of variables ti, . . . ,tq, . . .. Thus, we 
obtain a basis of the latter as follows: for a given Young diagram A, define the Schur function sx[t] 
by 

(1.13) SA[t] = (£|e^M|A;£) 

(by translational invariance it is in fact independent of i). In the language of Schur functions, the 
tq are (up to a conventional factor 1/q) the power sums, see section ri.2.31 below. 

We provide here various expressions of sx[t] using the free fermionic formalism. In fact, many of 
the methods used are equally applicable to the following more general quantity: 

(1.14) sx/^[t] = {fi;£\e^i*^\X;i) 

where A and ^ are two partitions. It is easy to see that in order for S;)^/^[t] to be non-zero, C A 
as Young diagrams; in this case Sx/^ is known as the skew Schur function associated to the skew 
Young diagram A/^. The latter is depicted as the complement of /i inside A. This is appropriate 
because skew Schur functions factorize in terms of the connected components of the skew Young 
diagram X/fi. 

Examples: = ti, sg = ^tf - t2, s^jj = ^tf + 12, s ^-p = ^tf - t^. 

2 2 5412 12 

■S jn = Si-i = ti, g j-jn = 24*1 + 2*1*2 + 2*2 ~ *1*3 " *4- 

1.2.2. Wick theorem and Jacobi-Trudi identity. First, we apply the Wick theorem. Consider as 
the definition of the time evolution of fermionic fields: 

In fact, (|1.12p gives us the "solution" of the equations of motion in terms of the generating series 
V(z), ^P*{z). 

Noting that the Hamiltonian is quadratic in the fields, we now state the Wick theorem: 
(1.16) {i\ ^P^M ■ ■ ■ V'.JOJV* M . . . ^l[t] \i) = det {i\ ^P^,ml[t] \i) 

^ l<p,q<n ■'^ 

Next, start from the expression ()1.14p of sa/^M- padding with zeroes A or ^ so that they have 
the same number of parts n, we can write 

sx/,[t] = {-n\ • • • ^^^_.e^Wv.t^„ . • • • rx^_^^^ \-n) 

and apply the Wick theorem to find: 

sx/jt]= det (-nlV,, „.ie-^Mv* \-n) 

It is easy to see that (— n| ipie^^^^ip* |— n) does not depend on n and thus only depends on j — i. Let 
us denote it 



(1.17) hk[t] = {l\e^i%l^^ |0) ^/ifc[i]z^ = (l|e^WV'*(z)|0) = e^'?>i*''^'' 

^ fc>o 

{k = j — i; note that = for k < 0). 



The final formula we obtain is 
(1-18) sx/f,[t] = det {hx^_^^_g+p[t]) 

or, for regular Schur functions, 

(1.19) sx[t]= det {hx^^,+p[t]) 

l<p,q<n 

This is known as the Jacobi-Trudi identity. 

By using "particle-hole duality" , we can find a dual form of this identity. We describe our states 
ten 



in terms of hole positions, parametrized by the lengths of the columns and yU^, according to 



Again the Wick theorem applies and expresses sx/^ in terms of the two point-function (m| ^/;*e^[*]^/^j \m) 

which only depends on i — j = k and is given by 

(1.20) 

ek[t] = (-1)'^ (-1| e^Wv_fc+i |0) J^efcitjz'^ = (-1| e^Wv(-z) |0) = e^9>i^"^)' 

The finally formula takes the form 
(1-21) Sa^W = ^<^et^^„ (cA^-M^^-g+pW) 

or, for regular Schur functions, 

(1.22) SA[t]= det (ex'^.g+p[t]) 

l<p,q<n \ ^ / 

This is the dual Jacobi-Trudi identity, also known as Von Nagelsbach-Kostka identity. 

1.2.3. Weyl formula. In the following sections ll.2.3ffT72.61 we shall fix an integer n and consider 
the following change of variable (this is essentially the Miwa transformation [79]) tq = ^ SILi -^i ■ 
The Schur function becomes a symmetric polynomial of these variables Xi, which we denote by 
sa(xi, . . . ,Xn), and we now derive a different (first quantized) formula for it. 

Due to obvious translational invariance of all the operators involved, we may as well set £ = n. 
Use the definition (|1.6|) of |A) and the commutation relations p.l2|) to rewrite the left hand side as 

{n\ e^W |A; n) = e^i>^ *i (n| rizi)riz2) ■ ■ ■ r{zn) |0) | 

where | means picking one term in a generating series. 

We can easily evaluate the remaining bra-ket to be: (we now use the £ = notation for the l.h.s.) 



(0| e^W I A) = e^^>i ^^=1 n {zi - Zj)\^r.+^,-^ 



l<i<j<n 



^2 ■■■^n 



Now write tq = \ X]j=i x'^j ^-^id note that e^^i>^ ^»=i ^» = nrj=i(-'- ~ ziXj) ^ . We recognize (part 
of) the Cauchy determinant: 



(0| e^W |A) = ^^^^<ij<n{'^ - ^i^i) ^ 



Y{i<j{Xi-Xj) 
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At this stage we can just expand separately each column of the matrix {1 — xiZj) ^ to pick the 
right power of Zj] we find: 

(1.23) (0|e^W|A) = ^%M^Ii^^^^ 

Remark 1: defined in terms of a fixed number n of variables, as in (|1.23p . s\{xi, . . . , Xn) has the 
following group-theoretic interpretation. The polynomial irreducible representations of GL{n) are 
known to be indexed by partitions. Then s\{xi, . . . , Xn) is the character of representation A evalu- 
ated at the diagonal matrix diag(xi, . . . , Xn)- Hence, the dimension of A as a GL{n) representation 
is given by sa(1^_^) = ni<i<i<n(^i " ^ " +i)/(j - i). 

n 

Remark 2: the more general Miwa transformation allows for coefficients: tq = ^ SiLi OjX?. In 
particular if we use minus signs, we get the notion of plethystic negation. Combining it with the 
usual negation of variables is equivalent to transposing Young diagrams: indeed it amounts to 
exchanging the and the e^-ft]. In other words, 

SAM = SA'[-e*] - etg := i-iy-Hq 

More generally, one defines the supersymmetric Schur function sx{xi, . . . , Xn/yi, ■ ■ ■ , ym) to be equal 
to sM where tq = \{Y.U < " TT=l{-y^Y)■ 



1.2.4. Schur functions and lattice fermions. Note that the change of variables = | Sj=i ■'^'j ^-Hows 
us to write 

1=1 q>l ^ 

So we can think of the "time evolution" as a series of discrete steps represented by commuting 
operators exp(j)^{xi). In the language of statistical mechanics, these are transfer matrices (and the 
existence of a one-parameter family of commuting transfer matrices exp (j)^ (x) is of course related 
to the integrability of the model) . We now show that they have a very simple meaning in terms of 
lattice fermions. 

Consider a two-dimensional square lattice, one direction being our space ^ + ^ and one direction 
being time. In what follows we shall reverse the arrow of time (that is, we shall consider that 
time fiows upwards on the pictures), which makes the discussion slightly easier since products of 
operators are read from left to right. The rule to go from one step to the next according to the 
evolution operator exp(/)_|_(x) can be formulated either in terms of particles or in terms of holes: 

• Each particle can go straight or hop to the right as long as it does not reach the (original) 
location of the next particle. Each step to the right is given a weight of x. 

• Each hole can only go straight or one step to the left as long as it does not bump into its 
neighbor. Each step to the left is given a weight of x. 

Obviously the second description is simpler. An example of a possible evolution of the system with 
given initial and final states is shown on Fig. H^a). 

The proof of these rules consists in computing explicitly (^| e^+^^^ |A) by applying say (|1.21|) for 
tq = ^x'', and noting that in this case, according to (jl.20|) . en[t] = for n > 1. This strongly 
constrains the possible transitions and produces the description above. 
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Figure 1. A lattice fermion configuration and the corresponding (skew) SSYT. 



1.2.5. Relation to Semi- Standard Young tableaux. A semi-standard Young tableau (SSYT) of shape 
A is a filling of the Young diagram of A with elements of some ordered alphabet, in such a way that 
rows are weakly increasing and columns are strictly increasing. 

We shall use here the alphabet {1,2, ... ,n}. For example with A = (5, 2, 1, 1) one possible SSYT 
with n > 5 is: 



1 


2 


4 


5 


5 


CO 


CO 
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5 











It is useful to think of Young tableaux as time-dependent Young diagrams where the number 
indicates the step at which a given box was created. Thus, with the same example, we get 







A 



So a Young tableau is nothing but a statistical configuration of our lattice fermions, where the 
initial state is the vacuum. Similarly, a skew SSYT is a filling of a skew Young diagram with the 
same rules; it corresponds to a statistical configuration of lattice fermions with arbitrary initial and 
final states. The correspondence is exemplified on Fig. [T]|^b). 

Each extra box corresponds to a step to the right for particles or to the left for holes. The initial 
and final states are and A, which is the case for Schur functions, cf (|1.13p . We conclude that the 
following formula holds: 



(1.24) 



sx{xi 



E n 

reSSYT{A,n) b box of T 



This is often taken as a definition of Schur functions. It is explicitly stable with respect to n in 
the sense that s\{xi, . . . , Xn, 0, . . . , 0) = sx{xi, . . . , Xn). It is however not obvious from it that sx 
is symmetric by permutation of its variables. This fact is a manifestation of the underlying free 
fermionic ( "integrable" ) behavior. Of course an identical formula holds for the more general case 
of skew Schur functions. 
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Figure 2. Underlying directed graphs for particles and holes. 



1.2.6. Non-Intersecting Lattice Paths and Lindstrdm-Gessel-Viennot formula. The rules of evolu- 
tion given in section 11.2.41 strongly suggest the following explicit description of the lattice fermion 
configurations. Consider the directed graphs of Fig. [2] (the graphs are in principle infinite to the left 
and right, but any given bra-ket evaluation only involves a finite number of particles and holes and 
therefore the graphs can be truncated to a finite part). Consider N on- Intersecting Lattice Paths 
(NILPs) on these graphs: they are paths with given starting points (at the bottom) and given 
ending points (at the top), which follow the edges of the graph respecting the orientation of the 
arrows, and which are not allowed to touch at any vertices. One can check that the trajectories of 
holes and particles following the rules described in section [1.2.4l are exactly the most general NILPs 
on these graphs. 

In this context, the Jacobi-Trudi identity ()1.19p becomes a consequence of the so-called Lind- 
strom-Gessel-Viennot formula [72l[36]. This formula expresses N{ii, . . . , ji, . . . , jn), the weighted 
sum of NILPs on a general directed acyclic graph from starting locations ii, . . . ,in to ending loca- 
tions j'l, ... ,in) where the weight of a path is the products of weights of the edges, as 

(1.25) N{ii,...,in]ji,...,jn) = detN{ip;jg) 

p,q 

More precisely, in Lindstrom's formula, sets of NILPs such that the path starting from i^ ends at 
jw{k) gst an extra sign which is that of the permutation w. This is nothing but the Wick theorem 
once again (but with fermions living on a general graph), and from this point of view is a simple 
exercise in Grassmannian Gaussian integrals. In the special case of a planar graph with appropriate 
starting points (no paths are possible between them) and ending points, only one permutation, say 
the identity up to relabelling, contributes. 

In order to use this formula, one only needs to compute N{i;j), the weighted sum of paths from 
i to j. Let us do so in our problem. 

In the case of particles (left graph), numbering the initial and final points from left to right, we 
find that the weighted sum of paths from i to j, where a weight Xi is given to each right move 
at time-step i, only depends on j — i; if we denote it by hj-i{xi, . . . ,Xn), we have the obvious 
generating series formula 

n ^ 

hkixi, . . . , Xn)z^ = TT 

■■^ 1 — ZXi 

k>0 i=l 

Note that this formula coincides with the alternate definition ()1.17p of hk[t] if we set as usual 
~ g Sr=i -^i • Thus, applying the LGV formula ()1.25p and choosing the correct initial and final 
points for Schur functions or skew Schur functions, we recover immediately (|1.18ll.l9p . 
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In the case of holes (right graph), numbering the initial and final points from right to left, we 
find once again that the weighted sum of paths from i to j, where a weight Xi is given to each left 
move at time-step i, only depends on j — i; if we denote it by ej-i{xi, . . . , Xn), we have the equally 
obvious generating series formula 

n 

Skixi, ... , Xn)z^ = JJ(1 + ZXi) 
k>0 i=l 

which coincides with (jl.20p . thus allowing us to recover (jl.21ll.22p . 

1.2.7. Relation to Standard Young Tableaux. A Standard Young Tableau (SYT) of shape A is a 
filling of the Young diagram of A with elements of some ordered alphabet, in such a way that 
both rows and columns are strictly increasing. There is no loss of generality in assuming that the 
alphabet is {1, . . . , n}, where n = |A| is the number of boxes of A. For example. 



1 


2 


6 


8 


9 


3 


4 






5 








7 









is a SYT of shape (5, 2, 1, 1). 

Standard Young Tableaux are connected to the representation theory of the symmetric group; the 
number of such tableaux with given shape A is the dimension of A as an irreducible representation of 
the symmetric group, which is up to a factor n! the evaluation of the Schur function s\ at tg = 5iq. 
Indeed, in this case one has H[t] = Ji, and there is only one term contributing to the bra-ket 
(A| e-f^[*l |0) in the expansion of the exponential: 

^A[5i.] = ^(A|jrio) 

In terms of lattice fermions, Ji has a direct interpretation as the transfer matrix for one particle 
hopping one step to the left. As the notion of SYT is invariant by transposition, particles and holes 
play a symmetric role so that the evolution can be summarized by either of the two rules: 

• Exactly one particle moves one step to the right in such a way that it does not bump into 
its neighbor; all the other particles go straight. 

• Exactly one hole moves one step to the left in such a way that it does not bump into its 
neighbor; all the other holes go straight. 

An example of such a configuration is given on Fig. O 

1.2.8. Cauchy formula. As an additional remark, consider the commutation of e^'^M and e^*t"l, 

where H*[u], the transpose of is obtained from it by replacing Jg with J_q. Using the 

Baker-Campbell-Hausdorff formula and the commutation relations (jl.lOp we find 

^H[t]^H*[u] ^ j:g>iqtgUg^H'[u]^H[t] 

or equivalents e'?^+(^)e'^-(?') = ^-l^e'^-(2^)e'^+(^) with (f>±{x) = Y,g>i ^J±q- 

If we now use the fact that the |A) form a basis of Tq, we obtain the Cauchy formula: 

(1.26) (0| e^We^'N |o) = ^ sx[t\sx[u\ = - x^y.Y^ = e^'/>i 



13 




1 — 

1 


ffl 


9 




6 

1 


B 




□ 


1 

) — 






1 


3 


2 


4 



(a) (b) 
Figure 3. A lattice fermion configuration and the corresponding SYT. 






(a) (b) 
Figure 4. (a) A plane partition of size 2x3x4. (b) The corresponding dimer configuration. 



with tn 



1 v^n 



q 



a i-ji=\ y 



1.3. Application: Plane Partition enumeration. Plane partitions are a well-known class of 
combinatorial objects. The name originates from the way they were first introduced |74j as two- 
dimensional generalizations of partitions; here we shall directly define plane partitions graphically. 
Their study has a long history in mathematics, with a renewal of interest in the eighties [9 
combinatorics, and more recently in mathematical physics [86]. 



m 



1.3.1. Definition. Intuitively, plane partitions are pilings of boxes (cubes) in the corner of a room, 
subject to the constraints of gravity. An example is given on Fig. Hl^a). Typically, we ask for the 
cubes to be contained inside a bigger box (parallelepiped) of given sizes. 

Alternatively, one can project the picture onto a two-dimensional plane (which is inevitably what 
we do when we draw the picture on paper) and the result is a tiling of a region of the plane by 
lozenges (rhombi with 60/120 degrees angles) of three possible orientations, as shown on the right 
of the figure. If the cubes are inside a parallelepiped of size a x b x c, then, possibly drawing the 
walls of the room as extra tiles, we obtain a lozenge tiling of a hexagon with sides a, b, c, which is 
the situation we consider now. 

Note that each lozenge is the union of two adjacent triangles which live on an underlying fixed 
triangular lattice. So this is a statistical model on a regular lattice. In fact, we can identify it with a 
model of dimers living on the dual lattice, that is the honeycomb lattice. Each lozenge corresponds 
to an occupied edge, see Fig. IHb). Dimer models have a long history of their own (most notably. 
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Kasteleyn's formula [5T] is the standard route to their exact solution, which we do not use here), 
which we cannot possibly review here. 



1.3.2. MacMahon formula. In order to display the free fermionic nature of plane partitions, we 
shall consider the following operation. In the 3D view, consider slices of the piling of boxes by 
hyperplanes parallel to two of the three axis and such that they are located half-way between 
successive rows of cubes. In the 2D view, this corresponds to selecting two orientations among the 
three orientations of the lozenges and building paths out of these. Fig. [5] shows on the left the 
result of such an operation: a set of lines going from one side to the opposite side of the hexagon. 
They are by definition non-intersecting and can only move in two directions. Inversely, any set of 
such NILPs produces a plane partition. 

At this stage one can apply the LGV formula. But there is no need since this is actually the case 
already considered in section 1.3.4. Compare Figs. [5] and [1} the trajectories of holes are exactly 
our paths (the trajectories of particles form another set of NILPs corresponding to another choice 
of two orientations of lozenges). If we attach a weight of Xi to each blue lozenge at step i, we find 
that the weighted enumeration of plane partitions in a a x 6 x c box is given by: 

Na,bA^l^ • • • ' ^a+b) = (0| \bx c) = Sbxcixi, Xa+b) 

where 6x c is the rectangular Young diagram with height b and width c. In particular the unweighted 
enumeration is the dimension of the Young diagram 6 x c as a GL{a + b) representation: 

i=ij=ik=i ■' 

which is the celebrated MacMahon formula. But the more general formula provides various refine- 
ments. For example, one can assign a weight of q to each cube in the 3D picture. It can be shown 
that this is achieved by setting Xi = q°-~^^~^ (up to a global power of q). This way we find the 
g-deformed formula 

a b c j+j+fc_i 

iVa,,c(g)=nnn i:k>.-2 

j=ij=ifc=i ^ 

Many more formulae can be obtained in this formalism. The reader may for example prove that 

Na,b,c= Yl sa(1,...,1)sa(1,...,1) 



A:Ai<c 



or that 



Na,b,c — det(l +TcxbTbxaTaxc) 

(where Tyxx is the matrix with y rows and x columns and entries (*), i = 0, — 1, j = 

0, . . . ,x — 1), as well as investigate their possible refinements, (for more formulae similar to the 
last one, see [3U]). Finally, one can take the limit o, 6, c — > oo, and by comparing the power of the 
factors 1 — in the numerator and the denominator, one finds another classical formula 

oo 

iVoo,oo,oo(g) = n(l-9")"" 
n=l 

Note that our description in terms of paths clearly breaks the threefold symmetry of the original 
hexagon. It strongly suggests that one should be able to introduce three series of parameters to 
provide an even more refined counting of plane partitions. With two sets of parameters, this is 
in fact known in the combinatorial literature and is related to so-called double Schur functions 
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Figure 5. NILPs corresponding to a plane partition. 



(these will reappear in section [5.2. 5p . The full three-parameter generalization is less well-known 
and appears in |lU7j . as will be recalled in section 14.3.21 

Remark: as the name suggests, plane partitions are higher dimensional versions of partitions, 
that is of Young diagrams. After all, each slice we have used to define our NILPs is also a Young 
diagram itself. However these Young diagrams should not be confused with the ones obtained from 
the NILPs by the correspondence of section [L2l 



1.3.3. Totally Symmetric Self- Complementary Plane Partitions. In the mathematical literature, 
many more complicated enumeration problems are addressed, see [97] • In particular, consider 
lozenge tilings of a hexagon of shape 2a x 2a x 2a. One notes that there is a group of transformations 
acting naturally on the set of configurations. We consider here the dihedral group of order 12 
which is consists of rotations of 7r/3 and reflections w.r.t. axis going through opposite corners of 
the hexagon or through middles of opposite edges. To each of its subgroups one can associate an 
enumeration problem. 

Here we discuss only the case of maximal symmetry, i.e. the enumeration of Plane Partitions with 
the dihedral symmetry. They are called in this case Totally Symmetric Self-Complementary Plane 
Partitions (TSSCPPs). The fundamental domain is a twelfth of the hexagon, see Fig. [6l Inside 
this fundamental domain, one can use the equivalence to NILPs by considering green and blue 
lozenges. However it is clear that the resulting NILPs are not of the same type as those considered 
before for general plane partitions, for two reasons: (i) the starting and ending points are not on 
parallel lines, and (ii) the endpoints are in fact free to lie anywhere on a vertical line. However the 
LGV formula still holds. For future purposes we provide an integral formula for the counting of 
TSSCPPs where a weight r is attached to every blue lozenge in the fundamental domain |28J. 

Let us call rj the location of the endpoint of the j*^ path, numbered from top to bottom starting 
at zero. We first apply the LGV formula to write the number of NILPs with given endpoints to be 

det{Ni^rj)i<i,j<n-i where Ni^r = (2j-r-i) ~ (1 + 'Tuy\y2i-r-i. Next we sum over them and 
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Figure 6. A TSSCPP and the associated NILP. 



obtain 

^n{T)= V det[(l+Tni)'<^]|™-i 0,-1 

0<ri<r2<---<r'n-i 

(1.28) „_i 

= n(l + ™.r E det«0|n_ 

j=l 0<ri<r2<-<r„_i 

We recognize the numerator of a Schur function; the summation is simply over all Young diagrams 
with n parts. At this stage we use a classical summation formula, s\{ui, . . . , Un^i) = nr=i^(-'- ~ 
Ui)~^ ni<i<j<n-i(l ~ 'UiUj)~^ , to conclude that 



n-1 . 2i — 1 



Uj -Ui TT (1 + TUi 



l<«<j<n— 1 «=1 



In 



n-1 . 2i-l 



where L-rn-i 2^-1 is now interpreted as picking the coefficient of a monomial in a powers series 
'lii^i 

around zero. 

This formula can be used to generate efficiently these polynomials by computer; in particular, 
we find the numbers 

A^„(l) = 1,2,7,42,429... 
which have only small prime factors. This allows to conjecture a simple product form: 

AT ( 1 ) = ff + = 

' J-J- (n + iV. n!(n + l)!...(2n-l)! 

i=0 

which was in fact proven in [2] . As a byproduct of what follows (sections 12.51 and 14. 4p , we shall 
obtain a (rather indirect) derivation of this evaluation. 

1.4. Classical integrability. The free fermionic Fock space is also important for the construction 
of solutions of classically integrable hierarchies. We cannot possibly describe these important ideas 
here, and refer the reader to [H] and references therein for details. Since an explicit example will 
appear in section 2, let us simply say a few general words. Recall the isomorphism $ 1— > {£\ e^'*] 
from J^£ to the space of polynomials in the variables tg (or equivalently to the space of symmetric 
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Figure 7. All TSSCPPs of size 1, 2, 3. 
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Figure 8. A configuration of the six- vertex model. 



functions if the tq are interpreted as power sums). The resulting function will be a tau- function of 
the Kadomtsev-Petiashvili (KP) hierarchy (as a function of the ig) for appropriately chosen |$). 
By appropriately chosen we mean the following. 

In the first quantized picture, the essential property of free fermions is the possibility to write 
their wave function as a Slater determinant; this amounts to considering states which are exterior 
products of one-particle states. Geometrically this is interpreted as saying that the state (defined 
up to multiplication by a scalar) really lives in a subspace of the full Hilbert space called a Grass- 
mannian. The equations defining this space (Pliicker relations) are quadratic; these equations are 
differential equations satisfied by {l\ e^^^ |$). They are Hirota's form of the equations defining the 
KP hierarchy. 

In section 2 we shall find ourselves in a slightly more elaborate setting, which results in the Toda 
lattice hierarchy. 



2. The six-vertex model 



The six vertex model is an important model of classical statistical mechanics in two dimensions, 
being the prototypical (vertex) integrable model. The ice model (infinite temperature limit of 
the six-vertex model) was solved by Lieb [69] in 1967 by means of Bethe Ansatz, followed by 
several generalizations |681 [70l [7T] . The solution of the most general six vertex model was given by 
Sutherland [99] in 1967. The bulk free energy was calculated in these papers for periodic boundary 
conditions (PBC). Here our main interest will be in a different kind of boundary conditions, the 
so-called Domain Wall Boundary Conditions. But first we provide a brief review of the six-vertex 
model. 
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Figure 9. Weights of the six- vertex model. 
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2.1. Definition. 

2.1.1. Configurations. The six-vertex model is defined on a (subset of the) square lattice by putting 
arrows (two possible directions) on each edge of the lattice, with the additional rule that at each 
vertex, there are as many incoming arrows as outgoing ones. See Fig. [8] for an example, and for 
two alternative descriptions: the "square ice" version in which arrows represent which oxygen atom 
(sitting at each lattice vertex) the hydrogen ions (living on the edges) are closer to, with the "ice 
rule" that exactly two hydrogen ions are close to each oxygen atom; and the "path" version in 
which one considers edges with right or up arrows as occupied, so that they form north-east going 
paths. Around a given vertex, there are only 6 configurations of edges which respect the arrow 
conservation rule, see Fig. [9l hence the name of the model. 



2.1.2. Weights. The weights are assigned to the six vertices, see Fig. [9l Thus the partition function 
is given by 

Z = (weight of the vertex) 

configurations vertex 

An additional remark is useful. With any fixed boundary conditions, one can show that the 
difference between the numbers of vertices of the two types c is constant (independent of the 
configuration). This means that only the product = C1C2 of their two weights matters. 

Let us denote similarly = aia2 and = 6162. One can write 



ai = ae 



02 = ae"-^""-^^ 



hi = be'^^^^y 



and consider that a, b, c are the weights of the vertices, while Ex, Ey are electric fields. In what 
follows, we shall consider by default the model without any electric field, where the Boltzmann 
weights are invariant by reversal of every arrow and ai = 02 = a, oi = 02 = a, ai = 02 = a; and 
sometimes comment on the generalization to non-zero fields. 

There is another way to formulate the partition function, using a transfer matrix. In order to set 
up a transfer matrix formalism, we first need to specify the boundary conditions. Let us consider 
doubly periodic boundary conditions in the two directions of the lattice, so that the model is defined 
on lattice of size M x L with the topology of a torus. Then one can write 



Z 



where Tl is the 2^ x 2^ transfer matrix which corresponds to a periodic strip of size L. Explicitly, 
the indices of the matrix are sequences of L up/down arrows. can itself be expressed as a 
product of matrices which encode the vertex weights; in the case of integrable models, we usually 
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denote this matrix by the letter R: 



(2.1) 



Then we have 
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(2.2) 



Tl — tro(-RoL • • • R02R01) 
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where means the matrix R actin; 



on the tensor product of i and j spaces, and is an 
additional auxiliary space encoding the horizontal edges, as on the picture (note that the trace is 
on the auxiliary space and graphically means that the horizontal line reconnects with itself). On 
the picture "time" flows upwards and to the right. 

The introduction of a vertical electric field amounts to multiplying the transfer matrix by an 
operator which commutes with it, of the form e^"^^ (S^ being the number of up arrows minus 
the number of down arrows). More interestingly, adding a horizontal field amounts to twisting 
the periodic transfer matrix: indeed, all the horizontal fields, using conservation of arrows at each 
vertex, can be moved to a single site, so that the transfer matrix becomes, up to conjugation by 



(2.3) Tl = troiRoL ■ ■ ■ Ro2Roi^) 

where the twist Q acts on the auxiliary space and is of the form 



2.2. Integrability. 



I X — q ^ 



X — X 



2.2.1. Properties of the R-matrix. Let us now introduce the following parametrization of the weights: 

a 

(2.4) b 

c = q-q~^ 

X, q are enough to parametrize them up to global scaling. Instead of q one often uses 



A 



a2 + 62 _ c2 



-1 



2ab 2 

In general, q or A are fixed whereas x is a variable parameter, called spectral parameter. It can be 
thought itself as a ratio of two spectral parameters attached to the lines crossing at the vertex. 

The matrix R{x) then satisfies the following remarkable identity: (Yang-Baxter equation) 



Rl2{x2/xi)Rl3{x3/xi)R23ix3/x2) = R23,{x3/ X2)Rl3{x3/ Xi)Ri2{x2/ Xi) 



This is formally the same equation that is satisfied by S matrices in an integrable field theory 
(field theory with factorized scattering, i.e. such that every S matrix is a product of two-body S 
matrices) . 
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Figure 10. Phase diagram of the six- vertex model. 



The i?-matrix also satisfies the unitarity equation: 



Ri2{x)R2i{x ^) = {qx-q ^){qx ^-q ^x) 




with X = X2/X1. The scalar function could of course be absorbed by appropriate normalization of 



2.2.2. Commuting transfer matrices. Consider now the transfer matrix as a function of the spectral 
parameter x, possibly with a twist: 



We thus have an infinite family of commuting operators. In practice, for a finite chain Tl{x) is a 
Laurent polynomial of x so there is a finite number of independent operators. 

Note that we could have used the more general inhomogeneous transfer matrix 



where now we have spectral parameters yi attached to each vertical line i and one more parameter 
xq attached to the auxiliary line. Then the same commutation relations hold for fixed yi and 
variable xq. 

As is well-known, the commutation of the transfer matrices is only one relation in the Yang- 
Baxter algebra generated by the so-called RTT relations. The latter lead to an exact solution of 
the model using Algebraic Bethe Ansatz |31j . 

2.3. Phase diagram. The phase diagram of the six-vertex model in the absence of electric field 
is discussed in great detail in chapter 8 of [1]. It can be deduced from the exact solution of the 
model using Bethe Ansatz after taking the thermodynamic limit. The physical properties of the 
system depend only on A = (g + g~^)/2, x playing the role of a lattice anisotropy parameter. We 
distinguish three phases, see Fig. (TO) 



R. 



(2.5) 



Tl{x) = tio{RoL{x) • • • Ro2{x)Roi{x)n) 




Tl{xo;xi, ...,xl)= tTo{RoL{yL/xo) • • • Ro2{y2/xo)Roi{yi/xo)Q) 
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Figure 11. Correspondence between (a) (A 
domino tilings. 



0) six-vertex (b) NILPs and (c) 



(1) A > 1: the ferroelectric phase. This phase is non-critical. Furthermore, there are no local 
degrees of freedom: the system is frozen in regions filled with one of the vertices of type 
a or b (i.e. all arrows aligned), and no local changes (that respect arrow conservation) are 
possible. 

(2) A < —1: the anti-ferroelectric phase. This phase is non-critical. This time there is a finite 
correlation length. The ground state of the transfer matrix corresponds to a state with zero 
polarization (in the limit A — > — oo, it is simply an alternation of up and down arrows). 

(3) —1 < A < 1: the disordered phase. This phase is critical. It possesses a continuum limit 
with conformal symmetry, and this limiting infra-red Conformal Field Theory is well-known: 
it is simply the c = 1 theory of a free boson on a circle with radius R given by 
A = — cos 7, < 7 < vr. 



2(1-7 A) ' 



The phase diagram in the presence of an electric field is more complicated, though the basic 
division into the three phases above remains. See |95[ for a description^ 



2.4. Free fermion point. Inside the disordered phase, there is a special point A = 0. We provide 
various representations of the six-vertex model which display the free fermionic behavior of this 
region of parameter space. 



2.4.1. NILP representation. It is tempting to try to interpret the "north-east going paths" of Fig. [8] 
as Non-Intersecting Lattice Paths. The problem is that they can touch at vertices. One way to fix 
it is to consider the slightly modified paths of Fig. [TTTb) The rule is to replace each vertex of (a 
with the corresponding dotted square of (b) and then patch together the latter to form the paths 
Note that the correspondence is no longer one-to-one: each vertex of type c\ corresponds to two 
possible local paths. 



"'^Note that the discussion of the phase diagram in |104) is incomplete. 

^Going from (a) to (b) amounts to combining the equivalences of [IT] and |105j . 
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The directed graph of the NILPs is the basic pattern 1^ repeated, with paths moving 
upwards and to the right, and with weights indicated on the edges. Comparing the weights we get 
the relations 

oi = a(55 (22 = 1 

6i = f3e 62 = 07 

Cl = (5 + e7 C2 = a/3 

Combining these we find that 0102 + 61^2 — ciC2 = 0, so the correspondence only makes sense at 
A = (and there are really only 4 parameters and not 5 as one might naively assume). 

2.4.2. Domino tilings. There is also a prescription to turn six-vertex configurations into domino 
tilings that is illustrated on Fig. [TlI c) |105j . As already mentioned, going from (b) to (c) is nothing 
but a slightly modified version of the bijection of [17] between NILPs and domino tilings. 

In order to understand the correspondence of Boltzmann weights, note that patching together 
the pictures of Fig. fTTTc) produces dominoes that span three dotted squares, for example 




In particular, one half of the domino is contained inside one square. This allows to classify dominoes 
into four kinds, depending on which half of the square it occupies (these are called north-, west-, 
south-, and east-going in [H]). Going back to Fig. [TlT c). we conclude that ai, 02, &i, 62 can be 
considered as the Boltzmann weights of the four kinds of dominoes. Furthermore, we have the 
relations 

Cl = 0102 + 6162 C2 = 1 
from which we derive as expected aia2 + 6162 — C1C2 = 0. 

Just as plane partitions are dimers on the honeycomb lattice, domino tilings can be considered 
equivalently as dimers on the square lattice. 

2.4.3. Free fermionic five-vertex model. The general five-vertex model is obtained by sending one of 
the a or 6 weights to zero while all other weights remain finite; in other words, one simply forbids 
one of the 6 types of vertices. For a discussion of the general five- vertex model , see for example |83j 
and in particular its appendix A. In the first part of this section, we choose to send both horizontal 
and vertical electric fields to minus infinity and a to zero, in such a way that oi becomes zero. In the 
representation in terms of north-east going paths, this amounts to forbidding crossings; however, 
these paths in general interact when they are close to each other. The paths become NILPs (i.e. 
they only interact through the Pauli principle) only if their weights are products over the edges, 
which implies that 6162 = C1C2. This leads us back to the model of section [2.4. 11 but with 5 sent to 
zero: what we get this way is the free fermionic five- vertex model, first discussed in [lOlj . 

If (5 = the NILPs of Fig. [Tllb) simply live on a regular square lattice, and of course at this 
stage we recognize the transfer matrix discussed in section 11.2.41 and illustrated on Fig. [1] (plain 
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Figure 12. From the five- vertex model to dimers or plane partitions. 




Figure 13. From the five- vertex model to dimers or plane partitions, dual version. 

lines). In section [L3l on plane partitions, it was also identified with the transfer matrix of lozenge 
tilings. To complete the circle of equivalences, we show on Fig. 1121 how to go from NILPs to either 
dimers on the honeycomb lattice or lozenge tilings, following Reshetikhin |93j . 

There is a second case which is worth mentioning (if only because it will reappear in section 
I5.2.5P : suppose instead that we send 62 to zero. This time the north-east going paths cannot go 
straight east any more. In this case it is natural to redraw all north-east moves with a right turn 
as straight lines (not just south-side-goes to east but also west-side-goes-to-north), and this way 
we recognize the dashed lines of Fig. [H with a slight modification: the whole picture is distorted in 
such a way that each path moves one step further to the right (so that north-west becomes north, 
and north becomes north-east). If we want these paths to be NILPs, we reproduce the weights of 
12.4.11 with 7 = 0. Finally the correspondence to lozenge tilings/dimers is illustrated on Fig. [T3l 

Note a difference between the models of lozenge tilings corresponding to these two versions of 
the free fermionic five-vertex models: the vertical spectral parameters fiow north-east in the first 
picture, whereas they flow north-west in the second picture. Ultimately, this is related to two 
possible inhomogeneous versions of Schur functions (double vs dual [double] Schur functions in the 
language of [5D]). See also the recent work |lU9j where these lozenge tilings are embedded in a more 
general square-triangle-rhombus tiling model. 

2.5. Domain Wall Boundary Conditions. Domain Wall Boundary Conditions (DWBC) were 
special boundary conditions which were originally introduced in order to study correlation functions 
of the six- vertex model [60]. However they are also interesting in their own right. 
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Figure 14. An example of configuration witli Domain Wall Boundary Conditions. 



2.5.1. Definition. DWBC are defined on a n x n square grid: all the external edges of the grid are 
fixed according to the rule that vertical ones are outgoing and horizontal ones are incoming. An 
example is given on Fig. [TH 

To each horizontal (resp. vertical) line one associates a spectral parameter Xj (resp. yj). The 
partition function is thus: 

n 

Z„(xi,...,x„;yi,....?/„) = ^ 

configurations i,j=l 

where w = a,b,c depending on the type of vertex (cf ()2.4p ). Here we do not allow any electric field 
for the simple reason that with DWBC (as with any fixed boundary conditions), using the same 
type of arguments as in the previous section, one can push the effect of the field to the boundary, 
where it only contributes a constant to the partition function. 

Remark: the (one-to-many) correspondence of section 12.4.21 sends DWBC six- vertex configura- 
tions to domino tilings of the Aztec diamond [56]. 



2.5.2. Korepin's recurrence relations. In [60], a way to compute Z„ inductively was proposed. It is 
based on the following properties: 



Zi=q-q 



-1 



Zn{xi, . . . , Xn] yi, ■ ■ ■ -yn) IS a Symmetric function of the {xj} and of the {yi} (separately). 
This is a consequence of repeated application of the Yang-Baxter equation (or equivalently 
of one of the components of the so-called RTT relations): 



(qyi+i/vi-Q ^yi/yi+i)Zn{- ■ ■ ,yi,yi+i, ■ ■ ■) = {qyi+i/yi-q ^yt/yi+i)^ 



25 



Xl 



1 


\ 1 


\ / 


\ / 




1 
























/ \ 


/ N 


1 


\ 






/ \ 


1 



yi 



Figure 15. Graphical proof of the recursion relation. 
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and similarly for the Xj. 

Zn multiplied by x^~^ (resp. is a polynomial of degree at most n — 1 in each variable 

xj (resp. yf). This is because (i) each variable say Xi appears only on row i (ii) a, b are 
linear combinations of x^^ , Xi and c is a constant and (iii) there is at least one vertex of 
type c on each row/column. 
The Zn obey the following recursion relation: 



(2.6) Zn{xi, . . . ,Xn]yi = Xl, . . . ,yn) 



{Q-Q ^)Y\.(l^^/^i ~ 1 JJ((??/j/2;i - '^Xi/yj)Zn-lix2,...,Xn;y2,---,yn 

i=2 



i=2 



Since yi = xi implies h[yi/xi) = 0, by inspection all configurations with non-zero weights 
are of the form shown on Fig. [TSJ This results in the identity. 



Note that by the symmetry property, Eq. (j2.6p fixes Zn at n distinct values of yi: Xj, i = 1, . . . , n. 
Since Z.„ is of degree n — 1 m y\, it is entirely determined by it. 
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2.5.3. Izergin's formula. Remarkably, there is a closed expression for Zn due to Izergin |40^ I39j . It 

is a determinant formula: 

(2.7) 

ni,j=iixj/yi-yi/xj)iQ^j/yi-Q ^Vi/xj) ^ ^ { q-q~^ 



Z„ = : , ; ^ det 



lli<i<j<ni^i/^j - Xi/xi){yi/yj - yj/yi) i,j=i-n \ {xj/yi - yi/xj){qxj/yi - q ^yi/xj) 

The hard part lies in finding the formula, but once it is found, it is a simple check to prove that it 
satisfies all the properties of the previous section. The symmetry under interchange of variables is 
evident from the structure of the formula, and the recurrence formula follows from looking at the 
zeroes outside the determinant and the poles inside the determinant: indeed, they must compensate 
each other for the result to be non-zero at say yi = xi, which immediately leads to expanding the 
determinant on the xi row and to the recurrence. 

2.5.4. Relation to classical integrability and random matrices. The Izergin determinant formula is 
curious because it involves a simple determinant, which reminds us of free fermionic models. And 
indeed it turns out that it can be written in terms of free fermions, or equivalently that it provides a 
solution to a hierarchy of classically integrable PDE, in the present case the two-dimensional Toda 
lattice hierarchy. We cannot go in any details here but provide a few remarks. 

Consider a function of two sets of n variables of the form 

(2.8) A(X)A(y) 
where if X = (xi, . . . ,a;„), Y = {yi, . . .,yn), then 

det(j){X,Y)= det (l){xi,yj) A{X) = FT {xi - xj) 

i,j=l,...,n 

l<i<j<n 

Certainly Izergin's formula ()2.7p is of this form (up to some prefactors and to — > x?, yj — > yf); 
but it is also the case of the Cauchy formula ()1.26p (under the form of the Cauchy determinant) 
and even of Weyl's formula ()1.23p (once divided by A(A); though usually one considers it as a 
function of the Xi only, which results in a tau- function of the KP hierarchy only). It also appears 
in problems of random matrices, and in particular in the Harish Chandra-Itzykson-Zuber integral 
[TUl [5H] (see |1U6| ). r„ is of course symmetric by permutation of variables in X, and in Y. 

We now make use of the following set of bilinear determinant identities: 

n+l m+l 

(2.9) det (j){X - Xi, Y) det ^{X' + Xi, Y') = ^ det (j){X' , Y' - y'j) det (j){X, Y + y'j) 
i=i j=i 

where X = {xi, . . . ,Xn+i), Y = (yi,...,y„), X' = {x[, . . . , x'^), Y' = {y[, . . . ,y'^^^). 

Using as in section 1 the Miwa transformation!^ ~ q X^ILi -^i^ ^g ~ q '^7=1 Ui^ ^^^'^ similarly for 
the primed variables, we find after standard contour integration tricks that (|2.9p can be rewritten 
in terms of the as 

'^'^ — t'^U,] .'lpE,>i(*,-t',K 
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n"^-"T„,[t - [u],s]Tn+i[t' + M, s']e^«>i 



— ^;"-™T„[t', s' - [v]]Tm+i [t, s + [v]]e^.>M-^^)-' 



2'Kiv 



"^Since we have only a finite number of variables, the correct prescription is to consider a symmetric polynomial 
in n variables as a linear combination of Schur functions with fewer than n rows. 
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where [u] = {u~^, u~^/2, . . . , /q, . . .), and where the integrals § du/{2Triu) have the meaning of 
picking the constant term of a Laurent series. Expanding this equation in powers of tq — t'^ and 
Sg — s'g results in an infinite set of partial differential equations satisfied by the Tn- They are the 
Hirota form of the two-dimensional Toda lattice hierarchy. Inside it, there are two copies of the KP 
hierarchy corresponding to varying only one set of variables (the tq or the Sq) and keeping n fixed. 
Tn is the tau-function of the hierarchy. 

In particular, if one expands to first order in ti — t[ and set m = n — 1, we find 
^oinA d d d d 

Oti OSi oil OSi 

which is a form of the Toda lattice equation. 

There is another representation which is particular useful for the homogeneous limit. Consider 
the Laplace (or Fourier, we are working at a formal level) transform of (p: 

<t>{x,y)= lldf,{a,b)e'^^ + yb 

Then one can write 

= lArr^A^V^ /•••/ flMa^A) . det (e^^«i) . .det (e^^^^i) 

n\A{X)A{Y) J J ^J^ ij=l,...,n jj=l,...,n 

This is formally identical to the partition function of a generalized two-matrix model with external 
fields for both matrices. 

Next, let us consider the homogeneous limit Tn{x, y) of such a function Tn{X, Y) where all Xj tend 
to X and all yi tend to y. Noting that detie'''''^ )/ A{X) ~ e^'^«"^A(a) where c„ = (Efc=i ^0"^ 
(by the usual trick of taking Xi = x + ie, e ^ 0) and similarly for the yi, 

Tn{x,y) = CnCn+i [ ■■■ / J] tf/i(a„ 6i)A(a)e^ ^"=1 «^ A(6)e^ ^i=i 

i=i 

This is a generalized two-matrix model with linear potentials. With an arbitrary potential, the 
partition function of such a model is known to be a tau-function of the two-dimensional Toda 
lattice hierarchy. In fact, we have the following fermionic representation, with notations similar to 
section [D 

Tnix,y) oc (n,n|e^'?>i'^"^+'i ^^-^-'i^ J (i^(a, (a)V'* (6))" |0, 0) 

Since we only have here linear potentials i.e. the primary times (the "ti"), we shall only recover 
the first equation of the hierarchy. Let us do so. First note the determinant formula 

rnix,y) = c2 . det // d^(a,6)aVe^« + y& = . det (^J^cl>{x,', 

ij=0,...,n-l J J ij=0,...,n-l \OX^ OyJ 

Of course the latter form could have been derived directly from (|2.8p . Next apply to either of these 
expressions the Desnanot-Jacobi determinant identity. More precisely, consider the matrix of size 
n + 1 and write that its determinant times the determinant of the sub-matrix of size n — 1 with last 
two rows and columns removed equals the difference of the two possible products of determinants 
of sub-matrices of size n with one row, one column, removed and the other row, other column 
removed (among the last two rows and columns). The result is: 

2 d d d d 

(2.11) n Tn+lTn-l = T„ — — r„ - ^T„— T„, 

OX ay ox oy 

(the factor of n? takes care of the c^). 
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Finally, in the special case that (^(x, y) only depends onx — y, then the previous formulae simplify. 
The measure o?//(a, b) is concentrated on a = — 6 and we have 

I\{X)I\{Y ) J J i,]=l,...,n ij=l,...,n 

which is a generalized one-matrix model with external field. In the homogeneous limit, Tn{x,y) 
becomes a function of a single variable t = x — y and we can write 

(2.12) ^n{t)= [ ■■■ [fl d^{a,)Aiafe* SLi «^ 

i=i 

which is a one-matrix model with linear potential. With an arbitrary potential, the partition 
function of such a model is known to be a tau-function of the Toda chain hierarchy. Writing as 
before 

rJi) = det / / du(a)a'+^e^^ = det ( -^-^(pit)] 

^ ^ (n!)2 i,j=0,...,n-l J J ' (n!)2 i,j=0,...,n-l \dt'+i^^ 7 

and applying the Jacobi-Desnanot identity we obtain 

(2.13) n Tn+lTn^l = TnT^ - T„ 

which is a form of the Toda chain equation. 



2.5.5. Thermodynamic limit. In [6T], the Toda chain equation ()2.13p above was used to derive the 
asymptotic behavior of the partition function of the six-vertex model with DWBC in two of its 
three phases: ferroelectric and disordered. These are the two phases where we expect the limit 
n — > oo to be smooth, as we shall explain. In this case, making the Ansatz that the free energy is 

2 r 

extensive, we plug the asymptotic expansion t„ ~ e~'^ ■' into the Toda chain equation and are left 
with a simple differential equation to solve: 

(2.14) /" = e'^f 

In fact, this is the simplest reasonable Ansatz that is compatible with (j2.13p . so that any solution 
of (j2.13p with a smooth large n limit will be governed by the differential equation (j2.14p . 

This is how the computation of the bulk free energy of the six-vertex model with DWBC for 
A > — 1 is performed in [61]. It is much simpler than the computation for PBC and the result is 
given in terms of elementary functions, and is thus different from that of PBC. Explicitly, we find 

{max(a, b) A > 1 
iab-^ |A|<1, g = > 0), x = e*(*+^) 



(2.15) zy™' 



cos -rvt/y 

Note the obvious interpretation of the result in the ferroelectric regime ~ a frozen phase where 
almost all arrows are aligned. 

In the anti- ferroelectric phase A < — 1 one expects a less smooth large n limit because the 
alternation of arrows that is favored in this phase will interact with the boundaries. This statement 
is made more precise in |105j , where matrix model techniques are used to compute the large n limit 
of ()2.12p (which essentially boil down to an appropriate saddle point analysis of it). And indeed 
one finds that logr„ has an oscillating term of order 1, which explains that the simple differential 
equation ()2.14p cannot account for the asymptotic behavior. 

In any case, in all phases one finds a result that is different from that of PBC. The explanation 
of this phenomenon is that the six-vertex model suffers from a strong dependence on boundary 
conditions due to the constraints imposed by arrow conservation. In particular there is no ther- 
modynamic limit in the usual sense (i.e. independently of boundary conditions). In |104j it was 
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Figure 16. From six- vertex to ASMs. 



suggested more precisely that the six-vertex model undergoes spatial phase separation, similarly to 
plane partitions [13] and other dimer models [53]. In other words, even far from the boundary of 
the system, the system loses any translational-invariance and the physical behavior around a given 
point is a function of the local polarization: as such, the model can have several (possibly, all) of 
the three phases coexisting in different regions. This was motivated by some numerical evidence, 
as well as by the exact result at the free fermion point A = 0, at which the arctic circle theorem 
|46j applies: the boundary between ferroelectric and disordered phases is known exactly to be an 
ellipse (a circle for a = b) tangent to the four sides of the square. So the apparent simplicity of the 
computation of the bulk free energy for DWBC conceals a complicated physical picture. 

Since then, there has been a considerable amount of work in this area. There has been more 
numerical work [T]. Some of the results of [UT] have been proven rigorously and extended using 
sophisticated machinery in the series of papers [315117] by Bleher et al. Finally, the curve separating 
phases has been studied in the work of Colomo and Pronko [HI [T2], and recently they proposed 
equations for this curve in the cases a = b, A = =bl/2 jl6j . The point A = 1/2 is of special interest: 
it corresponds to all weights equal, and is the original ice model. It is also the subject of the next 
section. 



2.5.6. Application: Alternating Sign Matrices. Alternating Sign Matrices are an important class of 
objects in modern combinatorics [77]. They are defined as follows. An Alternating Sign Matrix 
(ASM) is a square matrix made of Os, Is and -Is such that if one ignores Os, Is and -Is alternate 
on each row and column starting and ending with Is. For example, 

10 
10 
10-11 
10 

is an ASM of size 4. The enumeration of ASMs is a famous problem with a long history, see [U]. 
Here we simply note that ASMs are in fact in bijection with six- vertex model configurations with 
DWBC |62]. The correspondence is quite simple and is summarized on Fig. [TBI For example. 
Fig. [H] becomes the 4 x 4 ASM above. 

We can therefore reinterpret the partition function of the six-vertex model with DWBC as a 
weighted enumeration of ASMs. It it natural to set the weight of all zeroes to be equal (a = 6), 
which leaves us with only one parameter c/a, the weight of a ±1. In fact here we shall consider only 
the pure enumeration problem that is all weights equal. We thus compute A = 1/2 and q = e*'^/^, 
and then Xi = q, yj = 1 so that the three weights are w{xi/yj) = q — q~^. 

At this stage there are several options. Either one tries to evaluate directly the formula (j2.7p : 
since the determinant vanishes in the homogeneous limit where all the Xi or yj coincide, this is a 
somewhat involved computation and is the content of Kuperberg's paper |62j . 
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There is however a much easier way, discovered independently by Stroganov and Okada 
|85| . It consists in identifying Zn at q = e*'^/^ with a Schur function. Consider the partition 
= [n — l,n — l,re — 2,n — 2,...,1,1), that is the Young diagram 

n-l 

(2.16) A(") = 



s^(n){zi, . . . , Z2n) is a polynomial of degree at most n — 1 in each Zi (use (jl.24p ) and, satisfies the 
following 

2n 

(2.17) S^(„){zi,.. ■,Zj = q~'^Zi, ...,Zn)= W{Zk- Zi)Sy^(„-i){zi, . . . , Zi, . . . , Zj, . . .,Z2n) 



k=l 



where the hat means that these variables are skipped (start from ()1.13p . find all the zeroes as 
Zj = q^Zi and then set Zi = Zj = to find what is left). 

This looks similar to recursion relations ()2.6p . After appropriate identification one finds: 



Z„(xi, . . . , x„; 2/1, . . . , yn)l^,..,, = (-l)-("-i)/2(g - q-^r ^^y^)~^'"'^ 

i=l 

/ 2 2 2 2 2 2\ 

Sx^iQ a^i, . . . ,g x^,y^,. . . ,yj 

Note that Zn possesses at the point q = e*"^/^ an enhanced symmetry in the whole set of variables 
{qxi, . . . ,qxn,yi, ■ ■ ■ ,yn}- Finally, setting Xi = q~^ and yj = 1 and remembering that this will 
give a weig ht of {q - Q-i)" to each ASM, one concludes that the number of ASMs is given by 

■> (n) _ . _ , (n) . 

An = 3--("-^)/2.,(„) (1, . . . , 1) = 3-"("-i)/2 n - — '._7 



Simplifying the product results in 

n-l 

= TT 

(n + i) 



(2.18) ^.= ni^±^ = l,2,7,42,429... 



i=0 

which is a sequence of numbers we have encountered before! In fact, the first proof of formula 
(j2.18p . due to Zeilberger |1U2| . amounts to showing (non-bijectively) that the number of ASMs is 
the same as the number of TSSCPPs. 

These are the ASMs of size 1, 2, 3 (+ = 1, - = -1): 

+ 

+ + 

+ +0 

+ 00 +00 + + + 00+ 00 + 

+ 00+ +00 00+ +-+ +00 + 

00+ + 00+ +00 + + +00 
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(a) (b) 

Figure 17. Vertices of (a) the CPL model and (b) the FPL model. 



As a check, one can take n — > oo in ()2.18p . and using Stirling's formula one finds 

3^/3 



This is to be compared with (j2.15p for 7 = 2tt/3, where we find Zn"" 
agree considering Z„, = {\/^i/2)"' An- 



9i/8. The two formulae 



3. Loop models and Razumov-Stroganov conjecture 



3.1. Definition of loop models. Loop models are an important class of two-dimensional statis- 
tical lattice models. They display a broad range of critical range of critical phenomena, and in fact 
many classical models are equivalent to a loop model. The critical exponents, formulated in the 
language of loop models, often acquire a simple geometric meaning; and many methods have been 
used to study their continuum limit, including the Coulomb Gaz approach [82], Conformal Field 
Theory [5] and more recently the Stochastic Lowner Evolution \65\ \66\ [67] . Here we are of course 
more interested in their properties on a finite lattice (in relation to combinatorics) and in the use 
of integrable methods. 

We shall introduce two classes of loop models on the square lattice, which turn out to be both 
closely related to the six- vertex model. Then we shall discuss a very non-trivial connection between 
these two loop models (the Razumov-Stroganov conjecture). 

Let us first discuss common features of the two models. Their configurations consist of loops 
living on the edges of the square lattice. The most important feature is the non-local Boltzmann 
weight produced by assigning a fugacity of r to closed loops. Here r is a real parameter (usually 
called n, due to the connection with the 0(n) model - however for various reasons we avoid this 
notation here). This can be supplemented by possible local weights for the various configurations 
around a given vertex, see Fig. [T71 

3.1.1. Completely Packed Loops. Configurations of Completely Packed Loops (CPL) consist of non- 
intersecting loops occupying every edge of the square lattice, which produces two possibilities at 
each vertex, represented on Fig. [TTlfa). Besides the weight r of each closed loop, one can introduce 
a local weight of u for one of the two types of CPL vertices, say NE/SW loops. 

The model is known to be critical for |t| < 2, and its continuum limit is described by a theory 

2 

with central charge c = 1 — 6 ^^J__^^ , where r = 2 cos 7, < 7 < vr. 

3.1.2. Fully Packed Loops: FPL and FPL^ models. Configurations of Fully Packed Loops (FPL) 
consist of non-intersecting loops such that there is exactly one loop at each vertex, which results 
in the six possibilities described on Fig. [TWb). One can then give a weight of r to each closed loop 
(FPL model). However, one can do better: noting that the empty edges also form loops (dashed 
lines on the figure), one can put them on the same footing as occupied edges and assign them a 
fugacity too, say f . This more general model is usually called FPL^ model. 
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Figure 18. From six- vertex to FPLs. 



One reason that the FPL^ is interesting is the following: the FPL model is not integrable for 
T 7^ 1 (the very special case r = 1 is of interest to us and will be considered below). However, the 
FPL^ is integrable for r = f, that is if the two types of loops are given equal weights. This was 
shown using Coordinate Bethe Ansatz in [2T] and then rederived using Algebraic Bethe Ansatz in 

For generic values of r, f , the Coulomb Gaz approach provides non-rigorous arguments to identify 

2 

the continuum field theory, see [H], and allows to compute the central charge to be c = 3— 6 ^^J_^^ — 

~2 

6 ^^J_--j , where r = 2 cos 7 and f = 2 cos 7, < 7,7 < vr. In particular the FPL model has central 

2 

charge c = 2 — 6 ^^J_^^ , which is one more than the corresponding CPL model. 

We shall not discuss the possibility of adding local weights in detail. Let us simply note that 
even if we impose rotational invariance of local Boltzmann weights, we can introduce an energy 
cost for 90 degrees turns of the loops, which amounts to giving them a certain amount of bending 
rigidity. Such a model was studied numerically in [32] • 



3.2. Equivalence to the six-vertex model and Temperley Lieb algebra. 



3.2.1. From FPL to six-vertex. The relation between six-vertex model and FPL model is rather 
limited, so we treat it first. The limitation comes from the fact that one cannot assign an actual 
weight to the loops, so that we obtain a r = 1 model (with only local weights). The correspondence 
between configurations is one-to-one: starting from the six-vertex model side, one imposes that at 
every vertex, arrows pointing in the same direction should be in the same state (occupied or empty) 
on the FPL side. This forces us to distinguish odd and even sub-lattices, and leads to the rules of 

Fig. da 

For rotational invariance of the FPL weights one should have a = b. c/a then plays the role of 
rigidity parameter of the loops mentioned in section 13.1.21 

The rest of this section is devoted to the equivalence of CPL and six- vertex models. 



3.2.2. From CPL to six-vertex. An example is shown on Fig. [T9j 

Start from a CPL configuration. The (unoriented) loops carry a weight of r. A convenient way 
to make the latter weight local is to turn unoriented loops into oriented loops: each configuration is 
now expanded into 2* ^oops configurations with every possible orientation of the loops. The weight 
of a 90 degrees turn is chosen to be oj^^^^, where t = uj -\- uj~^. 
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Figure 19. From CPLs to six-vertex. 
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Finally we forget about the original loops, retaining only the arrows. We note that the arrow 
conservation is automatically satisfied around each vertex: we thus obtain one of the six vertex 
configurations. 




Note that if u = 1 all weights become rotationally invariant and a = 6, ci = C2. 

Finally, one checks that the formula A = — r/2 holds (equivalently q = — w), u playing the role 
of spectral parameter. In particular the critical phase |A| < 1 corresponds to |r| < 2. 

Remark: this construction only works in the plane. On the cylinder or on the torus we have a 
problem: there are non-contractible loops which according to the prescription above get a weight 
of 2. This issue will reappear in the section 13.2.41 under the form of the twist. It explains the 
discrepancy of central charges between 6- vertex model (c = 1) and CPL (c < 1). 



3.2.3. Link Patterns. In order to understand this equivalence at the level of transfer matrices, one 
needs to introduce an appropriate space of states for the CPL model. We now assume for simplicity 
that L is even, L = 2n. 

Define a link pattern of size 2n to be a non-crossing pairing in a disk of 2n points lying on 
the boundary of the disk. Strictly equivalently we can map the disk to the upper half-plane and 
"fiatten" link patterns to pairings inside the upper half- plane of points on its boundary (a line). We 
shall switch from one description to the other depending on what is more convenient. The points 
are labelled from 1 to 2n; in the half-plane, they are always ordered from left to right, whereas in 
the disk the location of 1 must be chosen, after which the labels increase counterclockwise. 

Denote the set of link patterns of size 2n by P2n- The number of such link patterns is Cn = ni^+iy. ' 
the so-called Catalan number. 
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Example: in size L = 6, there are 5 link patterns: 
5 5 5 








3.2.4. Periodic Boundary Conditions and twist. Suppose that we consider the CPL model with 
periodic boundary conditions in the horizontal direction, with a width of L = 2n. We can define a 
transfer matrix with indices living in the set of link patterns P2n as follows. Consider appending a 
row of the CPL model to a link pattern; this way one produces a new link pattern: 





The transfer matrix Tt^tt' is then the sum of weights of CPL rows such that the pattern vr' is turned 
into the pattern tt. The weights are calculated as follows: first one takes the product over each 
plaquette of the local weights; and then one multiplies by r to the power the number of loops that 
we have created. 

What is the precise correspondence between the space of link patterns C-Pjn and the space of 
spins which relates the transfer matrix of the six-vertex model and the newly defined one for 
the CPL model? 



We start from the equivalence described in the section [3.2.2i The basic idea is to orient the loops. 
So we start from a link pattern and add arrows to each "loop" (pairing of points). Forgetting about 
the original link pattern we obtain a collection of 2n up or down arrows, which form a state of the 
6-vertex model in the transfer matrix formalism. To assign weights it is convenient to think of the 
points as being on a straight line with the loops emerging perpendicularly: this way each loop can 
only acquire a weight of w^^/^, depending on whether it is moving to the right of to the left. For 
example, in size L = 2n = 4, 



1 2 3 4 U 1 2 3 



-j-l 234 -j-l 234 ~\~0J 

^ Ti T2 is + Ti [2 Ts i4 + ii T2 is T4 +^ ii [2 Ts 



1 2 3 



234 "^1234 1234 

]l [2 Ts [a + Tl [2 is T4 + il T2 Ts i4 +^ ^ ii T2 is T4 



There is only one problem with this correspondence: it is not obviously compatible with periodic 
boundary conditions. We would like to identify a loop from i to j, i < j and a loop from j to i + L, 
j < i + L. This is only possible if we assume that Ti+L= Tn [i+L= ii, i-e. we impose twisted 
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boundary conditions on the six- vertex model. In the notations of (j2.3|) the twist is = uj'^' : it 
corresponds to an imaginary electric field. 

This mapping from the space of link patterns (of dimension c„) to that of sequences of arrows 
(of dimension 2^" ) is injective; so that the space of link patterns is isomorphic to a certain subspace 
" . The claim, which we shall not prove in detail here but which is a natural consequence of the 
general formalism is that the transfer matrix (j2.3p of the six- vertex model with the twist defined 
above leaves invariant this subspace and, once restricted to it, is identical to the transfer matrix of 
our loop model up to this isomorphism, the correspondence of weights being the same as in section 
[3X2] (in particular, A = -t/2). 

The connection between CPL and six-vertex models is deep, in the sense that they are based 
on the same algebraic structure, the Temperley-Lieb algebra and the associated solution of the 
Yang-Baxter equation, but in different representations. This is what we briefly discuss now. These 
definitions will serve again when we study the quantum Knizhnik-Zamolodchikov equation (section 

m). 



3.2.5. Temperley-Lieb and Hecke algebras. The Temperley-Lieb algebra of size L and with param- 
eter r is given by generators Cj, i = l,...,L — 1, and relations: 

(3.1) = TCi ejCjiiej = Cj Cje-,- = ejCi \i - j\ > I 

It is a quotient of the Hecke algebra, i.e. the satisfy the less restrictive relations 

2 



(3.2) 



TCi 



j\ > 1 



Note that in Hecke (not Temperley-Lieb!), there is a symmetry Cj r — Cj. The Hecke algebra is 
itself a quotient of the braid group algebra: if r = —{q + q~^) {q is thus a free parameter), then the 
ti = q~^/'^ei + q^/'^ satisfy the relations 



titj 



tjti 



j| > 1 



We are interested in two representations of the Temperley-Lieb algebra. 

The first one is the representation on the "space of spins", that is the same space (C^)'^'^ of 
sequences of up/down arrows on which the six- vertex transfer matrix acts. It is given by making 
Cj act on the i^^ and {i + copies of in the tensor product, with matrix 



(3.3) 



/O 







-q 
1 







1 

-q~ 










It may be expressed in terms of Pauli matrices as 



(3.4) 



1 



= i<rfaf^i + a^a^, + A(af <t^i - 1) + h{at^, - af)) 



where h = ^{q — q ^) and the are the Pauli matrices at site i, 

The second representation of Temperley-Lieb can be defined purely graphically. We now assume 
as before that L is even, L = 2n. In order to define the action of Temperley-Lieb generators Ci on 
the space of link patterns Ci-*2n (vector space with canonical basis the | tt) indexed by link patterns) , 



it is simpler to view them graphically as Ci 



; then the relations of the Temperley-Lieb 
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algebra, as well as the representation on the space of link patterns, become natural on the picture; 
for example, we find 



ei 





62 



1 



1 



1 



6 



2 3 4 5 6 

As before, the role of the parameter r is that each time a closed loop is formed, it can be erased at 
the price of a multiplication by r. 

The two representations we have just defined are of course related: the transformation of section 
13.2.41 makes the representation on link patterns a sub-representation of the one on spins. 

Finally, we need to define affine versions of Temperley-Lieb and Hecke algebras. It is convenient 
to do so by starting from their non-affine counterparts and adding an extra generator p, as well as 
relations pcj = ej+ip, i = 1, . . . , L — 1, and = 1. Note that this allows to define a new element 
cl = pe-L-iP~^ = P~^^iP such that all defining relations of the algebra become true modulo L. In 
fact, a more standard approach would be to introduce only cl and not p itself. Adding p leads to a 
slightly extended affine Hecke/Temperley-Lieb algebra, which is more convenient for our purposes 
(see the discussion in [57]). 

In the link pattern representation, p simply rotates link patterns counterclockwise: 




In the spin representation, it rotates the factors of the tensor product forward one step and twists 
the last one (that moves to the first position) by fi. Once again, for these two representations to 
be equivalent, one needs to be of the form Q = uj°"' where cu = —q. 

As an application, consider the Hamiltonian Hi, obtained as the logarithmic derivative of the 
transfer matrix Tl{x) of ()2.5p evaluated at x = 1. We find that with periodic boundary conditions, 
it is simply given up to additive and multiplicative constants by 



Hl 



In the spin representation, using (j3.4p . we recognize in Hl the Hamiltonian of the XXZ spin chain 
(with the so-called ferromagnetic sign convention). So the Hamiltonian of the loop model, which 
has the same form, is equivalent to the XXZ spin chain Hamiltonian, but with twisted periodic 
boundary conditions, which in terms of Pauli matrices means that cr^_)_^ ~ , ,±2^± 



1 • 
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Figure 20. The CPL model on a cylinder. 



3.3. Some boundary observables for loop models. Here we consider the CPL model at r = 1 
with some specific boundary conditions which will play an important role since the observables 
we shall compute live at the boundary. Several geometries are possible and lead to interesting 
combinatorial results [1^, but here we only consider the case of a cylinder. 



3.3.1. Loop model on the cylinder. We consider the model of Completely Packed Loops (CPL) on a 
semi-infinite cylinder with a finite even number of sites L = 2n around the cylinder, see Fig. [20j It 
it convenient to draw the dual square lattice of that of the vertices, so that the cylinder is divided 
into plaquettes. Each plaquette can contain one of the two drawings and 

We furthermore set r = 1 (or q = e^'^*/^), that is we do not put any weights on the loops. There 
are no more non-local weights, and in fact plaquettes are independent from each other. So we 
can reformulate this model as a purely probabilistic model, in which one draws independently at 
random each plaquette, with say probability p for -y^^ and I — p for 

Finally, we define the observables we are interested in. We consider the connectivity of the 
boundary points, i.e. the endpoints of loops (which are in this case not loops but paths) lying 
on the the bottom circle. We encode them into link patterns (see section r3.2.3p . In the present 
context, they can be visualized as follows. Project the cylinder onto a disk in such a way that 
the boundaries coincide and the infinity is somewhere inside the disk. Remove all loops except the 
boundary paths. Up to deformation of these resulting paths, what one obtains is a link pattern. 
The probabilities of occurrence of the various link patterns can be encoded as one vector with c„ 
entries: 

where P2n is the set of link patterns of size 2n and is the probability of link pattern vr. 



3.3.2. Markov process on link patterns. We now show that can be reinterpreted as the steady 
state of a Markov process on link patterns. This is easily understood by considering a transfer 
matrix formulation of the model. As in section 13.2.41 let us introduce the transfer matrix: it 
corresponds here to creating one extra row to the semi-infinite cylinder, and encoding not the 
actual plaquettes but the effect of the new plaquettes on the connectivity of the endpoints. That 
is, Ttt^tt' (p) is the probability that starting from a configuration of the cylinder whose endpoints are 
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connected via the link pattern tt' and adding a row of plaquettes, one obtains a new configuration 
wliose endpoints are connected via the hnk pattern vr. This form a x matrix Ti{p). 

This transfer matrix is actually stochastic in the sense that 

(3.5) T^y{v) = l Vvr' 

which expresses the conservation of probability. This is of course a special feature of the transfer 
matrix at r = 1. Note that (|3.5p says that Tl{p)'^ has eigenvector (1, . . . ,1) with eigenvalue 1. 

The matrix Tl{p) has non-negative entries; it is easy to show that it is primitive (the entries of 
Tl{p)^ are positive). These are the hypotheses of the Perron-Frobenius theorem. Therefore, Tl{p) 
possesses a unique eigenvector with positive entries; the corresponding eigenvalue is positive 
and is larger in modulus than all other eigenvalues. Now the theorem also applies to Tl{p)^ and 
by uniqueness we conclude that the largest eigenvalue of Tl{p) and of Tl{p)'^ is 1. In conclusion, 
we find that the eigenvector with positive entries of Tl (p) , which with a bit of foresight we call ^ l 
again, satisfies 

(3.6) Tl(j))'^l = -^L 

(In fact the whole reasoning in the previous paragraph is completely general and applies to any 
Markov process, being up to normalization the steady state of the Markov process defined by 
Tl{p).) 

Two more observations are needed. Firstly, (|3.6p is clearly satisfied by the vector of probabilities 
that we defined in the previous paragraph (the semi-infinite cylinder being invariant by addition of 
one extra row); it is in fact defined uniquely up to normalization by (|3.6|) . This explains that we 
have used the same notation. 

Secondly, is in fact independent of p. The easiest way to see this is to note that p now 
plays the role of spectral parameter (explicitly, with the conventions of the six-vertex i?-matrix, 
p = g!^ri_g-i^ )- In particular we conclude that, as in the six- vertex model, \Tl{p),Ti{p')\ = 0, so 
that these transfer matrices have a common Perron-Frobenius eigenvector. 

3.3.3. Properties of the steady state: some empirical observations. We begin with an example in 
size L = 2n = 8. By brute force diagonalization of the stochastic matrix Tl one obtains the vector 
of probabilities: 
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We recognize some of our favorite numbers An, namely 7 and 42. 

In fact, Batchelor, de Gier and Nienhuis [3] conjectured the following properties for all system 
sizes L = 2n: 

(1) The smallest probability is and corresponds to link patterns with all parallel pairings. 

(2) All probabilities are integer multiples of the smallest probability. 

(3) The largest probability is An-i/A^, and correspond to link patterns which pair nearest 
neighbors. 

By now all these properties have been proven |23l I28| . as will be discussed in section HI 

3.3.4. The general conjecture. A question however remains: according to property 2 above, if one 
multiplies the probabilities by A^, we obtain a collection of integers. The smallest one is 1 and the 
largest one is An-i, but what can we say about the other ones? 

Recall that An also counts the number of six- vertex model configurations with DWBC. Further- 
more, we showed that there is a one-to-one correspondence between six- vertex model configurations 
and FPL configurations (cf section r3.2.ip . In this correspondence, the DWBC become alternating 
occupied and empty external edges for FPL configurations (in short, FPLs). Let us now draw 
exphcitly the 42 FPLs of size 4x4: 
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Interestingly, we find that the reformulation in terms of FPLs allows to introduce once again a 
notion of connectivity. Indeed, there are 2n occupied edges on the exterior square and they are 
paired by the FPL. We can therefore count separately FPLs with a given link pattern vr; let us 
denote this by A-^. The Razumov-Stroganov conjecture [UD] then states that 

A-n 

thus relating two different models of loops (CPL and FPL) with completely different boundary 
conditions. And even though both models are equivalent to the six-vertex model, the values of A 
are also different (they differ by a sign). 

The Razumov-Stroganov conjecture remains open, though some special cases have been proved, 
see below. 

The relation to the conjectured properties of the previous section is as follows. It is easy to show 
that if TT is a link pattern with all pairings parallel, then there exists a unique FPL configuration 
with connectivity vr. Thus the RS conjecture implies property (1). Furthermore, since all are 
integer, it obviously implies property (2). Property (3) however remains non-trivial, since even 
assuming the RS conjecture it amounts to saying that yl^r = ^n-i in the case of the two link 
patterns vr that pair nearest neighbors, which has not been proven. 



4. The quantum Knizhnik-Zamolodchikov equation 

We now introduce a new equation whose solution will roughly correspond to a double gener- 
alization of the ground state eigenvector of the loop model introduced above: (i) it contains 
inhomogeneities and (ii) it is a continuation of the original eigenvector to an arbitrary value of q, 
the original value being q = e^'^^^^. 

4.1. Basics. 

4.1.1. The qKZ system. Let L = 2n be an even positive integer. Introduce once again the -R-matrix, 
but this time rotated 45 degrees and which acts a bit differently than before. Namely, it acts on 
the vector space CP2n spanned by link patterns, in the following way: 

(4.1) i?.(.) = ^(z)i^:i^4±il^ ^ = l,...,L-l 

q — q 

where we refer to section [3.2.51 for the definition of e^. (t){z) is a scalar function that we do not need 
to specify explicitly here, see section 14.1.21 for a discussion. Redrawing slightly the latter as = 
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>, and similarly the identity as 1 




we recognize the two (rotated) CPL plaquettes. 



Note that in this section, it is convenient to use spectral parameters z that are the squares of our old 
spectral parameters x. Indeed, using the equivalence to the six-vertex model described in section 
13. 2. 4^ which amounts to the representation ()3.3p for the Temperley-Lieb generators (acting on the 
i*^ and (i + 1)*^* spaces), we essentially recover the i2-matrix of the six- vertex model, cf (|2.1|2.4p . 
after the change of variables z 
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Vx'^^'^R{x)x where V permutes 



provided one performs the following transformation: R{z 
the factors of the tensor product, and k = diag(0, 1,-1,0). 

The Yang-Baxter equation can be rewritten as follows: 

Ri{z)Ri+i{zw)Ri{w) = Ri+i{w)Ri{z w)Ri{z) i = I, . . . ,L - 1 

We also require that (l){z)(j){l/ z) = 1, so that Ri satisfies the unitary equation: 

Ri{z)Ri{l/z) = 1 

Consider now the following system of equations for ^'i,, a vector- valued function of the zi 
{i = l,...,L-l) 

(4.2) Ri{zi+i/zi)'iiLizi, . ..,zl) = ^'l(2i, . . .,Zi+i,Zi, ... ,zl) 

(4.3) P"^^'l(2;i, . ..,Zl) = K^Il{z2, . . .,ZL,SZi) 

where it is recalled that p rotates link patterns counterclockwise!! and k is a constant that is needed 
for homogeneity, s is a parameter of the equation: if one sets s = q^^^~^^^ with k = 2 (technically, 

this is the dual Coxeter number of the underlying quantum group [7^(51(2))), then i is called the level 
of the qKZ equation. (|4.2|) is the exchange equation, which is ubiquitous in integrable models. (|4.3|) 
has something to do with our (periodic) boundary conditions, or equivalently with an appropriate 
affinization of the underlying algebra. 

This system of equations first appeared in |96j in the study of form factors in integrable models. 
It is not what is usually called the quantum Knizhnik-Zamolodchikov {qKZ) equation; the latter 
was introduced in the seminal q-deformation of the Knizhnik-Zamolodchikov (KZ) 

equation (gKZ is to quantum affine algebras what KZ is to affine algebras). The usual qKZ equation 
involves the operators Si, which Si can be defined pictorially as 

Si = 



-1. 



where the empty box is just the "face" graphical representation of the i?-matrix (dual to the 
"vertex" representation we have used for the six- vertex model): 



+ 



1 



and the spectral parameters z to be used in Si are as follows: for the box numbered j, zj/zi if j > i 
or Zj/{szi) if j < i. Loosely, Si is the "scattering matrix for the i*^ particle". 



^We use p ^ in the l.h.s. because (p ^'I'L)7r ~ 9 



p(t) 
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Alternatively, 5^ can be expressed as a product of Ri and p: 

Si = Ri-i{Zi-i/{sZi)) ■ ■ ■ R-2{z2/{sZi)Ri{zi/{sZi))pRL^l{zL/zi) ■ ■ ■ Ri+i{Zi+2/ Zi)Ri{zi+i/ Zi) 

The quantum Knizhnik-Zamolodchikov equation can then be written 

(4.5) Si{zi, . . . . . . , Zi, . . . , zl) = . . . , s Zi, . . . , zl) 

It is a simple exercise to check using ()4.4p that the system ()4.2H4.3p implies the qKZ equation 
()4.5p . Naively, the converse is untrue. However one can show that, up to some linear recombinations, 
a complete set of (meromorphic) solutions of ()4.5p can always be reduced to a complete set of 
solutions of (IMOp . see [H]- 

4.1.2. Normalization of the R-matrix. In the usual setting of the qKZ equation, the normalization 
4>{z) of the i?-matrix is chosen such that it satisfies an additional constraint, the crossing symmetry, 
as well as certain analyticity requirements. There is however a certain freedom in choosing this 
normalization. Indeed, consider a solution ^ of the system (|4.2|4.3p and redefine 

^= n /(^j-M)^ 

^<i<j<L 

where f{z) is some scalar function. Then ^ satisfies the same system of equations: the first 
equation remains formally identical, but with the normalization of the i?-matrix modified from 
(p{z) to (j){z) = <p{z) f {z) / f {1 / z) . For the second equation to remain true imposes a constraint on 
the function f{z), so that we find 

f{z)/f{l/z) = 4>{^)/Hz) 
f{sz) = f{l/z) 

With reasonable analyticity conditions on 0/(/>, this system of equations can be solved (e.g. if s < 1, 
writing i{z)/^{z) = fo{z)/foil/z), one finds /(z) = (11^=0 /o(sV^) 0^=1 /o(5'^))~')- 

In what follows we shall only be interested in homogeneous polynomial solutions of the system 
()4.2H4.3p . In this case the most convenient normalization is to set </> = !, which we shall adopt from 
now on. 

Note that iterating (|4.3p L times implies that k^s'^ = 1 where d is the degree of ^. 

4.1.3. Relation to affine Hecke algebra. The role of the affine Hecke algebra (and even double affine 
Hecke algebra) was emphasized in the work of Cherednik [12] in relation to Macdonald polynomials 
and the quantum Knizhnik-Zamolodchikov equation, and in the work of Pasquier [87] on integrable 
models. In particular, its use in the context of the Razumov-Stroganov conjecture was advocated 
by Pasquier [89] . 

Start from the qKZ system (j4.2H4.3p . and rewrite it in such a way that the action on the finite- 
dimensional part (on the space of link patterns) is separated from the action on the variables (on 
the space of polynomials of L variables). Start with (j4.2p : after simple rearrangements it becomes 

(4.6) (r - ei)^L = -{qzi - q-^Zi+l)^^^L 

where di = j-^^{si — 1) and Si is the operator that switches variables Zi and so that the 
l.h.s. only acts on the polynomial part of ^ whereas the r.h.s. only acts on link patterns. 

The operators r — Cj := — {qzi — q~^Zi+i)di acting on the space of polynomials form a representa- 
tion of the Hecke algebra (with parameter r); i.e. they satisfy the relations of ()3.2p . Equivalently, 
the ej satisfy them and form a second set of generators. 
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Figure 21. From a Young diagram to a link pattern: in this example, from the 
partition (2,1,1) to the pairings (1,2), (3,8), (4,5), (6,7). 

As to (|4.3|) . it is already in a separated form; and one can add an extra operator on the space 
of polynomials, the one that appears in the r.h.s. of (j4.3|) . that is the cyclic shift p of spectral 
parameters zi Z2 ■ ■ ■ zl s zi (with the multiplication by k included). The Cj together 
with p generate a representation of the affine Hecke algebra. We can write formally: 

(4.7) a^L = h'^L 

(4.8) p-^^L = P^L 

We can now interpret (j4.7l4.8p as follows: we have on the one hand a representation of the affine 
Hecke algebra on CP2n.) the space of link patterns (with generators Cj and p); and on the other 
hand a representation of the same algebra on C[zi, . . . the polynomials of L variables (the 
Cj and p). provides a bridge between these two representations: it is essentially an invariant 
object in the tensor product of the two, that is it provides a sub-representation of the space of 
polynomials (explicitly, the span of the ^jr) which is isomorphic to the dual of the representation 
on the link patterns. By dual we mean, defined by composing the anti-isomorphism that sends the 
Cj to themselves and p to p"^ (but reverses the order of the products) and transposition. The role 
of the anti-isomorphism is most evident when we try to compose equations (|4.7l4.8p and find that 
iiCj^ = CjCi^ (in fact, it is often advocated to make one of the sets of operators act on the right, 
though we find the notation too cumbersome to use here). 

So the search for polynomial solutions of (j4.2|4.3p is equivalent to finding certain irreducible 
sub-representations of the action of affine Hecke on the space of polynomials. 

Remark: the direct relation between qKZ and representations of an appropriate affine algebra 
only works for the An series of algebras i.e. affine Hecke. For more complicated situations such 
as the Birman-Wenzel-Murakami (BWM) algebra, it does not work so well because one cannot 
separate the two different actions [88J. 

4.2. Construction of the solution. On general grounds, we only expect polynomial solution for 
integer values of the level. Here we shall only need a solution at level 1, that is s = q^. 

4.2.1. qKZ as a triangular system. We shall build this solution in several steps. First, we use a 
"nice" property of our basis of link patterns, that is, the fact that (|4.6p can be written as a triangular 
linear system in the components of ^l. This requires to define an order on link patterns, which 
is most conveniently described as follows. Draw once again link patterns as pairings of points on 
a line and consider the operation described on Fig. [2TJ It gives a bijection between link patterns 
of size 2n and Young diagrams inside the staircase diagram 1" = (n — 1, n — 2, . . . , 1). Then the 
order is inclusion of Young diagrams. The smallest element, corresponding to the empty Young 
diagram, is denoted by 0"; it connects i and 2n + 1 — i (note that it is one of the link patterns 
with all pairings parallel, which correspond to the smallest probability 1/A„ in the loop model). 
Consider now the exchange equation (j4.6p and write it in components; we find two possibilities: 
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• i and i + 1 are not paired. Then we find that ()4.6p only involves ^'tt, and implies that 
qzi — divides ^'tt, and furthermore ^'tt/ {q Zi — q~^ZiJ^i) is symmetric in the exchange 
of Zi and Zj+i. 

• i and i + 1 are paired. Then 

[q Zi - q~^ Zi+i)di'^ = ^ 

that is it involves the sum over preimages of vr by viewed as acting on the set of link 
patterns. It turns out there are two types of preimages of a given vr: in terms of Young 
diagrams, there is the Young diagram obtained from vr by adding one box at i, i + 1 (which 
is always possible unless vr is the largest element); and there are other Young diagrams that 
are included in vr. So we can write the equation 

(4.9) ^'TT+one box at = {q Zi - q~^ Zi+i)di^ ^ - ^ ^ ^> 

it' ClT 

ei-7r'=7r 

which has the desired triangular structure and allows to build the ^'tt one by one by adding 
boxes to the corresponding Young diagram. However there is no equation for ^'q" . In fact 
this triangular system can be explicitly solved |54j (see also [19]) in the sense that every ^'^^ 
can be written as a series of operators acting on \I'on . We shall not use this explicit solution 
here. 

From the discussion above, we find that all we need is to fix ^'o". We use the following simple 
observation, which generalizes the first case in the dichotomy above: 

If there are no pairings between points i, i+ 1, . . . , j in vr, then Y\i<p<cq<j 

(qzp-q ^Zq) 

divides ^tt. 
(the proof is by induction on j — i)- 

We now make the "minimality" assumption that in the case of 0", these factors form all of ^'o", 
so that we find 

(4.10) ^on= Yl ilZi-q-^j) Yl {(l~^Zj-qZi) 

'^<'i<j<n n+l<i<j<2n 

where we recall that the system size is L = 2n. 

If a solution of the qKZ system with such a ^'o" exists, then it will be the only solution of degree 
n(n — 1) (up to multiplication by a constant). It remains however a non-trivial fact that with such 
a choice of ^'o", the construction above is consistent (i.e. independent of the order in which one 
adds boxes), and that ()4.3p is satisfied, with s = q^ . This is the subject of the next section. 

Example: we find in size L = 6 

^ _ {q^zi - Z2) {q^zi - Z3) {(f'Z2 - 2:3) {q^Zi - Z5) {(f'Zi - zq) {(f'z^ - zq) 



5 {q^Zi - Z2) {q'^Z3 - Zj) {q'^Z5 - Zq) 8,8 6 6 

W6,^^;X4= 7 [ziZ2Z3q + ziZ2Z4q - ziz^z^q - Z2Z3Z4q 

A ^ q^ 



ZlZ2Z5q'^ - ZxZ2Zio(f' + Z3Z!^Z^(^ + Z^Z^Z^f^ + Z^Z^Z^C^ + 22^5^60'^ " ^3^5^6 " Zi^Z'^Zi^) 



[c^Z2 - Z3) [(^Z2 - Zi) {q^zs - Zi) {q^zi - Z5) {q^zi - zq) {q^z^ - zq) 



q' 
45 




[q^zi - Z2) {q'^zs - Z4) {q^zs - Z5) {q^z^ - Z5) {q'^Zi - Zq) {q^Z2 - zq) 



{q^Z2 - zz) {q^Zi - Z5) {q^zi - zq) 



{ziZ2Z3q 
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Z2Z3Z4q'^ - Z2Z3Z5q^ - ZiZ^Z^q*^ 



ziZ2ZQq'^ - ziZsZQq*^ + Z2Z4Zriq^ + z^z^z^q'^ + Z2Z3ZQq'^ + ziz^zeq"^ + ziz^ZQq'^ - Zi^z^z^) 



4.2.2. Consistency and Jucys-Murphy elements. Let us go back to the qKZ system ()4.7|4.8p . In 
order to show that this system of equations is consistent, one must simply simply check that one 
obtains the same representation of affine Hecke on both sides (up to duality). Thus, we shall 
now check that the center has the same value in both representations. This will fix k and more 
importantly s. 

Let us first define a new set of elements of the affine Hecke algebra (Jucys-Murphy elements), 
namely 

yi = ti--- tL-ip~^tl^ ■ ■ ■ t~}^ i = l,...,L 

For this section we shall use "knot-theoretic" drawings to explain the equations we write. It is easy 
to see that if the ti = q~^/'^ei + q^/"^ are crossings in the knot theoretic sense (this definition being 
itself the skein relation for the Jones polynomial), then the yi have an equally simple description: 

where the second picture is embedded in a strip that is periodic in the horizontal direction. Then 
it becomes graphically clear that the yi commute: 



ViVj = 



One can also show that the symmetric polynomials of the y[ = q~''yi form the center of the 
affine Hecke algebraH Thus, assuming that we have an irreducible representation such that we can 
diagonalize the action of the y-, the (unordered) set of their L eigenvalues is independent of the 
eigenvector and characterizes the representation. 

Let us now apply this to our two representations. First, let us make the yi act on link patterns. 
They have an obvious common eigenvector: the link pattern 1" which pairs neighbors 2i — 1 and 





^This extra factor of comes from our normalization of the ti. In the Hecke algebra they have no preferred 
normalization, and one usually chooses it to make this factor disappear, but in Temperley-Lieb there is one (which 
makes the skein relation rotationally invariant), that we have chosen here, and it is unfortunately different. 
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2i, i = 1, . . . ,n. Indeed, we find: 



yd' 




-3/2 



7 D 



2 3 



t 

i.e. that the action of yi is equivalent to a Reidemeister move I (with negative orientation for i 
even and positive orientation for i odd). A move I multiphes the Jones polynomial by —q~^^^ (this 
amounts to the calculation tiCi = —q~'^^'^ei), so that we obtain 

(4.11) = -q^-^y^/h^ 

On the polynomial side, define similarly = ■ ■ ■ i^^pii-i ■ ■ - U, where ii = q~^/'^ei + q'^^"^ . 
One can show that there is an order on monomials (see [89]) such that the operators yi are upper 
triangular. Thus it suffices to evaluate their diagonal elements in the basis on monomials. The 
action of is most easily computed on monomials Hi^i 'h^ with non-increasing powers Aj+i < A^. 
In this case we find 

L L 

(4.12) iji JJ z^' = -«;s^»g3("+3-*) JJ + lower terms 

i=l i=l 

Comparing the two expressions (j4.1ip and (j4.12p . we find that s = q^ and A = (Aj) is, up to a 
global shift, equal to A^") = (n — 1, n — 1, . . . , 1, 1, 0, 0). This is the minimal degree solution that 
we shall consider - the degree being n{n — 1), as assumed earlier. We also deduce k = 

The latter computation looks very similar to a computation in nonsymmetric Macdonald poly- 
nomials [73]. This is no coincidence. In fact the eigenvectors of iji are exactly nonsymmetric 

Macdonald polynomials, see [50]. If we restrict ourselves to leading terms j^^^ z^^^'^^ which are in 
bijection with link patterns by labelling closings/openings in decreasing order, e.g. 



TT= I — * — i — i — \ — i — I — \ 1 — > A(7r) = (3,2,3, 1,2, 1,0,0) 
12345678 ' ' V'^'''^'/ 

then their span is exactly the span of the ^'7^, even though they are distinct from the ^'7^ - there is 
a triangular change of basis, so that only is itself a nonsymmetric Macdonald polynomial. 



4.2.3. Wheel condition. In [89], a particularly attractive characterization of the span of the ^'7^(2:1, . . . , 
71" S -fbn; is obtained. It will give an independent check of some of the results of the last section, as 
we shall see. 

Consider the following subspace of homogeneous polynomials: 

(4.13) M„ = {/(zi, . . . , Z2n), deg,^ / < n - 1 : /(. . . , z, . . . , q^z, . . . , gS, . . .) = 0} 

that is polynomials that vanish as soon as an ordered triplet of variables form the sequence 1, g^, g^. 
This vanishing condition is a so-called wheel condition, see [l9] for more general ones. We have 
also imposed a bound of n — 1 on the degree in each variable, which is a little bit more restrictive 
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than what is really needed: imposing on the (total) degree deg / < n{n — 1) would suffice, but this 
condition is less convenient for our purposes^ 

It is easy to check that ^'o" £ more interestingly, the action of r — Cj = —{qzi — q~^Zi+i)di 
preserves M„ (check the vanishing property by discussing separately the three cases: (i) neither i 
or i + 1 are in the triplet: trivial; (ii) both i and i + 1 are in the triplet: then the prefactor vanishes; 
and (iii) one of them is: then use the property for the original triplet and the one with i replaced 
with i + 1). Since the can be built from ^'o" by action of the r — Cj, one concludes that they 
are all in Af^Q In fact, note that the action of p also leaves Mn invariant on condition that s = q^, 
so that Mn is a representation of the affine Hecke algebra. 

We have an inclusion of the span of the ^'yr inside M„. In order to prove the equality, we shall 
state certain useful properties. 

4.2.4. Recurrence relation and specializations. First, we need the following simple observation. Ac- 
cording to the dichotomy of section 14. 2. H if vr is a link pattern such that the points i and z + 1 
are unconnected, then qzi — q~^Zi+i divides ^'7^, so that ^^(z^+i = q'^Zi) = 0. But if i and i + 1 
are connected (we shall call such a pairing of neighbors a "little arch"), then we can use the wheel 
condition to say that ^7r(-2^i+i = q'^Zi) vanishes when zj = q^Zi, j > i + 1, or zj = q~'^Zi, j < i, so 
that it is of the form n}=i(-2i ~ l'^ ^j) Ylj=i+ii^i ~ where by degree consideration <I> can 

only depend on the zj with j ^ i,i + 1. It is already clear that ^ is in M„,_i and so is a linear 
combination of the entries of ^l-~2', but there is better. 

Call ipi the mapping from Pl~2 to Pl which inserts an extra little arch at (z,i + 1). Then we 
claim that 

vr Im ifi 

J=l jti+2 ^ = '^^^ 

"^■jr'izi,. . . ,Zi-.i,Zi+2,- ■ ■ ,Zl) 

A direct proof of this formula is particularly tedious. Let us instead use the following trick. Consider 
the following cyclic invariance property: the formulae ()4.14p for different values ofi = l,...,L — 1 
can be deduced from each other by application of ()4.3p . In fact, one can even extend this way 
(|4.14p to i = L; being a little careful of the factors of s = that appear here and there, we obtain: 

2n-l 

(4.15) "^^U^'M = q'^ZL, ...,zl) = Yl - q^Zj)^^,{z2, . . . , zl-i) 

i=2 

where ipi adds an arch between L and 1 (on top of all other arches in the half-plane picture of link 
patterns). So, to prove (j4.14p . it is enough to prove (j4.15p . 

We shall use the construction of section 14.2.11 Note that the mapping is particularly natural 
via the bijection to Young diagrams, since it becomes the embedding of the set of diagrams inside 
1"-^ to that of dia grams inside 1". In particular 0" ^ is sent to 0"; so the first check is compatibility 
with ()4.10p . i.e. that when vr = 0" in ()4.15p we obtain in the r.h.s. the correct prefactors times ^'gn-i . 
The second check is that the formula ()4.9p which allows to compute the components of ^l, when 



(4.14) ^^{...,Zi+i=q^z, 



^In fact, one can show that having a partial degree greater than n — 1 while retaining a total degree of n{n — 1) 
would lead to a contradiction with our representation of affine Hecke using (|4.f 2p . 

'''a more visual proof of this result is given in |f08] . where it is noted that since + 1) and (i + f , i + 2) cannot 
be both pairings, ^L{zi,Zi+i = q^Zi+i, Zi+2 = q*Zi) = and then use (|4.2[l to permute the arguments and get the 
general wheel condition. 
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restricted to diagrams inside produces the components of This is rather obvious 

graphically: the action of Cj is compatible with the natural inclusion of diagrams (it can be defined 
directly in terms of Young diagrams). Furthermore, in order to remain inside 1""^, one must never 
add boxes at (1, 2) or (L — 1, L). So we never affect the variables zi, zl and up to a shift of indices 
i ^ i — 1, this is just the procedure in size L — 2. This is enough to characterize entirely the r.h.s. 

of Km\ . 

Next, we need two more closely related properties: 

• If / G M„, then 

(4.16) fiq''^^- ■ ■ , q~^'^") = V(ei) Dyck paths increments ^ / = 

(a Dyck path increment is a sequence {ei)i=i^,,,^2n of ±1 such that Yl'i=i = Si=i ^« — 
for all j < 2n). So a polynomial in M„ is entirely determined by its values at these Cn 
specializations. 

• We have the specializations: 
(4.17) 

(_l)n{n-l)/2(g _ ^-l)n{n-l)^|7r| ^. ^ gign(7r(i) - i) , i = 1, . . . , 2n 

otherwise 



(|7r| is the number of boxes in the correspondence with Young diagrams of Fig. l21|) . In other 
words, these values are the entries in the basis of the ^t^. 

The first point is shown explicitly in appendix C of [M! and we shall not reproduce the proof here; 
the second one is simply an application of the recurrence formula ()4.14p , by removing little arches 
one by one (see also [28]). 



4.2.5. Wheel condition continued. Using the properties of the previous section, we can now conclude 
that dimM„ < Cn and that the are independent polynomials, so that M„ is the span of the ^t^. 

One could use more representation-theoretic arguments to prove the equality of M„ and of the 
span of the ^'tt, for example based on the fact that as a representation of the affine Hecke 
algebra, is irreducible for generic q (in fact, it corresponds to the rectangular Young diagram with 
2 rows and n columns), but this is not our philosophy here. 

The wheel condition has the following physical interpretation (borrowed once again from |89j). 
Consider the case q = ±1. Then the qKZ equation ()4.5p becomes the ordinary KZ equation at level 
1. It is most conveniently expressed in the spin basis: 

OZi , Zi Zj 

where cTjj- exchanges spins i and j. The KZ equation is now in its usual form for correlation 
functions in s[(2)-^, except for a trivial change of ciij — 1/2 to Oij + 1, which is the same as 
^ — > ^nj<j(^i ~ ^j)^^^ (which is needed since we want a polynomial solution and not the usual 
one involving square roots). The solution is a product of two Vandermonde determinants: 

(4.18) = (-l)*^^<^-"'=-'"^=+^A(z,)„,=+A(^,)a.=- 

(ai,...,aL) G {+,-}^ 
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This is simply the ground state wave function of two species of free fermions ^l'±, i.e. with notations 
similar to section 1 

K^t;a,M, Z2n) = (0, 0| ^aAzi) ■ ■ ■ i^aAZn) \n, n) 

In this language the wheel condition is just the Pauli exclusion principle (if three fermions sit at 
the same spot, at least two of them must be of the same kind!). Of course the g = ±1 case is trivial 
in the sense that starting from say ^o"j one can apply the exchange equation (|4.2|) which simply 
says in this limit, using Ri = — cjj^j+i, that exchanging variables is the same as exchanging spins 
(up to a fermionic minus sign), which immediately results in ()4.18p . 

For a general value of q, one has a q-deformed version of this system of free fermions in which 
the vanishing condition of the wave function occurs not at coinciding points but when ratios of 
coordinates are ordered successive powers of q^. The usual fermionic statistic is turned into a 
non-trivial braid group statistic. 

4.3. Connection to the loop model. In general, the two problems of diagonalizing the transfer 
matrix of the CPL loop model and finding solutions of the gKZ system are unrelated. However 
there is exactly one value of q where a solution of gKZ does in fact provide an eigenvector of the 
transfer matrix (with periodic boundary conditions). This is when the parameter s = 1, which 
here occurs when q = e^'^^/^ (other sixth roots of unity are possible but they are either trivial or 
give the same result as the one we picked). In this case note that ()4.3p becomes a simple rotational 
invariance condition. Furthermore the real gKZ equation ()4.5p becomes an eigenvector equation 
for the scattering matrices: 

Si{zi, zl)^'l(zi, ...,zl) = "ifiizi, ... ,zl) 

These scattering matrices do not involve any extra shifts of the spectral parameters, and as is 
well-known in Bethe Ansatz, are just specializations of the inhomogeneous transfer matrix. Indeed 
if we define Tl{z; zi, . . . , zl) to be simply 

Tl = z \ I [ •••| 

Zl Z2 ZL 

(with periodic boundary conditions), where the box represents the i?-matrix evaluated at the ratio 
of vertical and horizontal spectral parameters, then observe that Si{zi., . . . , zl) = T^izi; zi, . . . , zl). 
By a Lagrange interpolation argument, we conclude that 

Tl(z;zi, . . . ,zl)^'l(zi, ...,zl) = *l(2;i, ■ ■ ■ ,2;l) 

i.e. ^' 2, is up to normalization the steady state of the inhomogeneous Markov process defined by 
r^(z; Zl, . . . , Zl). In order to recover the original homogeneous Markov process, one simply sets all 
Zi = 1. 

Note that the reasoning above is valid for any solution of gKZ, not just the one that was discussed 
in section 14.2.11 But there is no point in considering higher degree polynomial solutions since at 
q = e^'^*/'^ they will simply be the minimal degree of section Pl.2.1l times a symmetric polynomial. 

We now show how to apply the qKZ technology to prove some of the statements formulated in 
section [3l 

4.3.1. Proof of the sum rule. As a simple application of the above formalism, we explain how 
to recover the "sum rule" at q = e^'^*/^. Let us explain what we mean by that. We start by 
noting that the normalization of ^ is fixed once the value of ^'o" is specified; in particular, if all 
Zi = 1, using (j4.10p we find ^'o" = 3"("~i)/2^ to be compared with the (conjectured) probability 
1/An associated to 0" , cf section [3.3.31 In other words, with this normalization, one should have 
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^neP2n = 3"(""^)/2A„. This is what we are going to show now. In fact, following |23j, we are 
going to show more generally the inhomogeneous sum rule that X^^gp^^ ^7r(-zii • • • , Z2n) is equal 
to the Schur function s^(„){zi, . . . , Z2n), where A^") is defined in ()2.16p . that is up to a constant 
the partition function of the six- vertex model with DWBC at q = e*'^/^. Instead of the inductive 
method of [23], i.e. the comparison of the recurrence relations of sections 14.2.41 and 12.5.21 we shall 
here proceed directly. 

Define the covector v with entries = 1, tt £ P2n- The stochasticity property of a matrix is 
the fact that v is left eigenvector with eigenvalue 1; and it is satisfied at r = 1 by all the and 
henceforth by Ri{z): 

vei = v =^ vRi{z) =v i = l,...,L-l, q = 

Applying this identity to (|4.2|) . we immediately conclude that v ■'^ = '^^^p^^ is a symmetric 
polynomial of its arguments zi, . . . , z^. 

STeP2„ • • • ' -^2n) is a symmetric polynomial of degree n(n — 1) which satisfies the wheel 

condition of (j4.13p . We claim that this defines it uniquely up to normalization (see |108| for a 
similar claim in a more general fused model). The simplest proof is to use once again property 
(j4.16p i.e. that polynomials in M„ are characterized by the specializations ((7"*^^, . . . j*?"*^^")- But 
for a symmetric polynomial, all these specializations reduce to only one. This proves the claim. 

Next, one notes that s-^{n){z\^ . . . , Z2n) also satisfies these properties. The degree bound is clear 
from the definition ()1.24p of Schur functions; as to the wheel condition, it follows from formula 
()1.23p by noting that if Zk = q^Zj = q'^Zi, the sub- matrix of the matrix of the numerator with 
columns {i,j, k} is of rank 2 and thus the determinant of the whole matrix vanishes. 

Finally, to fix the normalization constant, note that thanks to (j4.17p . all components of ^ 
vanish at {zi,...,Z2n) = (q'^ , ■ ■ ■ ,q~^ ,q, ■ ■ ■ ,q) except ^'o"(g"^ . . . , g, . . . , g) = S"-^""^)/^; on 
the other hand, using ()2.17p . we find the same value s^(n) {q~^, • • • , q~^, q, ■ ■ ■ ,q) = 3"-("-i)/2, Wg 
conclude that 

^7r(^l, • • ■,Z2n) = S^(n){zi, . . . , Z2n) 

Remark: from the representation-theoretic point of view, v is an invariant element (under affine 
Hecke action) in (CP2n)*; and X^^gp^^ ^vr is its image in C[zi, . . . , zl]. This shows that at q = e^'^*/^, 
the representation of affine Hecke on CP2n (resp. the span of the ^tt) is not irreducible (though it 
remains indecomposable), since it has a codimension one (resp. dimension one) stable subspace. 

4.3.2. Case of few little arches. Note a remarkable property of recurrence relation ()4.14p : it can 
only decrease the number of little arches! So, if one considers the subset of link patterns with a 
given maximum number of little arches, that subset is closed under these relations. 

The link patterns that possess only two little arches (which is the minimum) are the link pattern 
0" and its rotations, i.e. the n link patterns for which all the pairings are parallel and for which we 
know explicitly the components: they are given by the action of p (cyclic rotation of variables plus 
shift of s) on ^'o" which is (|4.10p . 

These are the smallest components in the homogeneous limit. They correspond via the Razumov- 
Stroganov conjecture to a single FPL. As such they are clearly the "easiest" components. 

It is natural to look at what happens when one considers more little arches. Link patterns 
with three little arches are of the form of Fig. [22j In fact, the case of three little arches was first 
investigated in [29] on the other side of the Razumov-Stroganov conjecture: there, the problem of 
the enumeration of FPLs with such a connectivity was solved, with the remarkable result that the 
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Figure 22. Link patterns with 3 arches. 



number of FPLs with hnk pattern (a, b, c) is simply equal to the number of plane partitions inside 
an a X 5 X c box, that is given by the MacMahon formula (|1.27p . The proof is bijective and we 
shall not reproduce it here. Let us just comment on it. The crux of the bijection is the observation, 
known as "de Gier's lemma" , that in an FPL with given connectivity many edges are fixed by the 
simple requirement of the connectivity of the endpoints. Removing all these fixed edges we obtain 
a simpler enumeration problem, typically of plane partitions. Now the little arches play a key role 
in the sense that they are the ones that limit application of de Gier's lemma. In other words, the 
more little arches there are, the fewer fixed edges in FPLs, and therefore the more complicated the 
enumeration is. Once again, the number of little arches is a measure of complexity. 

Closing this philosophical parenthesis, let us go back to the gKZ equation and try to compute 
the entries ^a,b,c of corresponding to link patterns with three series of a, b, c nested arches. 
This is performed in |lU7j for q = e^'^'^l'^ and we generalize it here to arbitrary q. 

Up to rotation and use of (|4.3p . one can always assume that one of the little arches, say the 
one that is part of the series of a nested arches, is between L and 1. Let us further rename the L 
variables as follows: they become (g~^7i, . . . , q~'^'ya+b, oti, . . . , ab+c, 'f'Pi, • • • , q^f^a+c), see Fig. [22l 
Then according to the results of section 14. 2. H 

^aAc = n {QO(i-q~^aj) ]J {q Pi - Q'-^ Pj) Yi ((l^i - 

l<i<j<b+c l<i<j<a+c l<i<j<a+b 

b{a+b-l)+c{a+3b-l)+c^ /_ . x (^+6)(^+6-l)+c(t.+(c-l) 
q V / ^a^b,c 

where ^a,b,c is a polynomial that is symmetric in the {oi}, in the and in the {7j}. The powers 
of q and sign have been chosen carefully so that when we rewrite recurrence relation ()4.14p for 
Zi = ab+c, -2^1+1 = in terms of ^a,b,c, everything cancels out and we are left with 

a+b 

(4-19) ^a,b,c\Pi=Cb+c = n("^+'= ~ 7fc)^aAc-l 

k=l 

where the arguments (3i, Ob+c ^-re missing from ^a,b,c-i- Now since ^a,b,c is symmetric in the a^, 
this relation provides ^a,b,c at 6 + c values of (3i; and it is not hard to check that ^a,b,c is of degree b 
in Pi . Thus it is entirely determined by these values (and can for example be obtained by Lagrange 
interpolation) . 

To complete the (implicit) calculation of ^a,b,c one must provide an initial condition. At c = 0, 
the link pattern (a, b, 0) is nothing but the base link pattern 0", and in this case we find 

a b 

(4.20) ^aAO = nn(":'-^*) 

i=ij=i 

Some explicit formulae for $a,fe,c are given in |107j . In fact they are "triple Schur functions" in the 
sense of multi-Schur functions of |63j . 
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■ ai 

(a) (b) 
Figure 23. The three sets of spectral parameters associated to lozenge tilings. 



From the discussion above, it is natural to try to put non-trivial Boltzmann weights on lozenge 
tilings of a hexagon of sides a, b, c in such a way that their partition function inside an a x b x c 
box coincides with the inhomogeneous component $a,fe,c- Such a model is introduced in |lU7j . We 
point it out here because we believe this model is of some interest and might deserve further study. 
The spectral parameters live on the medial lattice, which is the Kagome lattice, see Fig. [23]^a). 
Each lozenge has a weight which is the difference of spectral parameters crossing at the center of 
the lozenge, with the sign convention a — /J,/? — 7, a — 70 

One can show that the integrability of the underlying model on the Kagome lattice implies that 
the partition function is a symmetric function of the of the {/3j} and of the {71}. We choose 

not to do so here and refer instead to the appendix of |107j for a "manual" proof of symmetry. 

Finally, it is easy to see that the partition function of lozenge tilings with such weights satisfies 
the recurrence relation (j4.19p . as illustrated on Fig. [23]^b); and that when c = 0, there is a unique 
tiling of the a x b parallelogram, resulting in (|4.20p . Note the similarity with the properties of 
the partition function of the six-vertex model with DWBC (symmetry and recurrence) found by 
Korepin. 

Thus, <J>a,6,c and this partition function coincide. Finally, if = 1, /3j = 7^ = (f', using the 
fact that the number of lozenges of each orientation is fixed and equal to ah, 6c, ca, we find: 



a,b,c 



a,b,c 



= (1 

-iQ 
-- (f 
(<?■ 



-2\ab 



2\ca 



N, 



a,b,c 



^—l\^ab+bc+ca ^a{c-b) ^_-Yy 



(a+c) (a+c- 1)+ (a+6) {a+b-l) 

(t.+ c)(6+e-l) . (c+a)(e+a-l) . (g + i)) (a + t- 1) 



''Na,b,c 

a+b~l)+c{a+'ib-l)+c'^ ^_ 



p27ri/3 



(a + i))(a + i)-l) + c(i)+(c-l) 
T) 2 



1 / -1 \ n(n — 1) 



'(-1)- 



a,b,c 



,27ri/3 



We conclude that ^'a,6,c/^0" = Na,b,c, as expected. 

The case of four little arches is similar to that of three arches and is sketched in |1U7| . On the one 
hand, the recurrence relations are once again enough to determine all ^-j^ uniquely. On the other 

^This is only a convention because the total number of lozenges of each orientation is conserved. The choice that 
we made does not respect the Z3 symmetry of the model but turns out to be convenient. 
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hand, the enumeration of FPL configurations with four Httle arches is doable (after appropriate 
use of the so-called Wieland rotation |100j ) and can be once again reduced to a lozenge tiling 
enumeration [30]. So the same strategy applies. Starting with five little arches, however, it fails on 
both sides: the recurrence relations are not sufficient to determine the corresponding components 
^'tt, and the enumeration of FPLs becomes non-trivial due to an insufficient number of fixed edges. 

4.4. Integral formulae. We now want to show that, using the formalism of the qKZ equation 
allows to prove the properties discussed in 13.3.31 as well as to reconnect the three models that 
we have found in which the same numbers An appear. A particular useful tool to exploit these 
solutions of qKZ has been introduced in [281 192j: it consists in writing integral formulae for them. 

4.4.1. Integral formulae in the spin basis. There are various types of known integral formulae for 
solutions of the quantum Knizhnik-Zamolodchikov equation. Here we are only interested in a very 
specific one, which only exists at level 1 and is obtained by (g-)bosonization, see |45j . 

The following formula is valid in even size L = 2n: 

(4.21) ^x%„_,(^i,...,^L) = (g-<?-^r n ii^-^'^^j) 

^<i<j<L 

n"=0 ni<i<a,("^^ - Ua,<^<LiQWi " Q'' 

where the contours surround the Zj counterclockwise, but not the q~'^Zi. The reader must be warned 
that we have changed the labelling of the entries of ^: the index, instead of being a sequence of 
L spins, is the ordered sequence of locations of plus spins. In other words the correspondence is 
(ai, . . . ^ {ao < • • • < «n.-i} = {i : ai = +}. 
For the case of odd size see [92] . 

4.4.2. The partial change of basis. Next we would like to write similar expressions for the entries of 
^ in the basis of link patterns. Unfortunately this would require to invert the change of basis from 
link patterns to spins that was described in section I3.2.4t and there is no simple explicit formula 
for the inverse (this is the famous problem of the computation of Kazhdan-Lusztig polynomials 
|52| : in this case there is a combinatorial formula [M] but it is not obvious how to combine it with 
integral formulae). 

Instead we shall do here something more modest: we introduce an intermediate basis between 
spins and link patterns. We skip here the details, which can be found in [28] (see in particular 
the appendices) and [27]. The bottom line is the introduction of a slight improvement of formula 
(lOB : 

(4.22) *ao,...,a„_i(^l,...,^L) 



(note the subtle modification 
the denominator). Here (oq. 




Jl {qzi-q ^Zj) 

l<i<j<L 

: the factors of disappeared and an inequality became strict in 
. ..,a„_i) is an arbitrary non-decreasing sequence of integers. The 
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connection with the entries of ^ in the bases of link patterns and spins is as follows: 

*ao,...,a„_i / V 'i/ ^ ao—eo,...,a„^i—e„-i 

£0.---,en-l6{0,l} 

*ao,...,a„_i = [ n ^#{^:i<«£<i}-0-i+l)/2 j ^TT 

7rGP2n \ i,i paired in tt / 

where Uk is the Chebyshev polynomial of the second kind evaluated at — r: Uk = - — q-'g-i • 

Note an interesting property of the second change of basis, which is that the coefficients are 
polynomials in r with integer coefficients. 

A second related property appears when one takes the homogeneous limit = 1 in the integral 
formula (j4.22p . It is convenient to send the two poles (1 and (/~^) to zero and infinity respectively 
by the homographic transformation U£ = {w£ — l)/{qw£ — q~^), so that: 

C""' • • • ' ^ / " / n 27rz n {Um-Ui){l+TUm+U(,Um) 

where r = —q — q~^, and the contours surround zero. This can be rewritten 



(1,...,!)= Y\_ {Um- Ue){l+TUm + UeUm)\ 



T71-1 



Wn'i 'ili=o "i 

" 0<i<m<n~l 

As shown in [27], one can find a subset of sequences {aQ, . . . , Un-i) such that the matrix of change 
of basis to link patterns is upper triangular with I's on the diagonal. Combining the two properties 
above, we reach the important conclusion that, normalized so that ^'o" = 1, the entries of the 
level 1 polynomial solution of qKZ in the homogeneous limit Zi = 1 are polynomials of r with 
integer coefficients. In fact, these coefficients are non-negative, but no combinatorial meaning is 
known for them. We shall however see in the next section that we can interpret coefficients of 
certain partial sums of entries in terms of plane partitions. 

Example: in size L = 6, we find the entries of ^ to be, with the same ordering of link patterns 
as in sections 13.2.31 and 14. 2. H the polynomials 1, 2r, r^, r^, r + r^. 



4.4.3. Sum rule and largest component. The integral formulae above allow to compute various linear 
combinations of the ^t^. Here, we concentrate on a single quantity of particular interest (which is 
considered in [28l [20]). Consider the linear combination Fn{t) = X^^gpj^ Cvr(t)^'7r(l) • • • > 1), with 
the normalization = !> ^-^d where the coefficient C7r(t) is most simply expressed in terms of the 
corresponding Young diagram by assigning a weight of t or t~^ to the boxes of even/odd rows in 
the complement inside 1", see Fig. [231 Three special cases are of particular interest: F{1) is just 
the sum X^^gp^^^ of all link patterns; F{0) is just ^'in, one of the two largest components; and 
F{t)\^n-i, the leading term of F(t), is just ^p(in-), the other largest component. 

We then claim that '^^^p^^ CTr{t)'^.^ can be naturally expressed in the intermediate basis as 
Z^eo e„-ie{o 1} '^'^i-eo,3-eiv,2n+i-e„_i (this is a simple calculation of the corresponding en- 
tries of the matrix of change of basis), which results in an integral formula. In the homogeneous 
limit we obtain 

n-l 

Fnit) = Y\_{'^+tUi) Y\. - Ul){l + TUm + UlUm)\YY:'-l 

i=0 0<e<m<n-l 
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Figure 24. Weight of a component expressed in terms of the associated Young diagram. 



In fact one can simply set = to get 

n.-l 

(4.23) i^„(t) = + (1 +rii£) {Um-Ui){l+TUm + UiUm)\YYi-i^2i-i 

i=l l<e<m<n~l 

which looks strikingly similar to the formula for the weighted enumeration TSSCPPs ()1.29p . 

There is however an important difference. Whereas (jl.29p just contains a product of functions 
of one variable times an antisymmetric function of the Ui (which, ultimately, comes from the free 
fermionic nature of the model), ()4.23p does not possess any particular symmetry w.r.t. exchange of 
its variables, due to the factors 1 + rUm + ugUm (which come from the g-deformed Vandermonde 
product in the original integral formula). One can however antisymmetrize formula ()4.23p . as 
conjectured in [2^ and then subsequently proved in |103j and in a slightly stronger version in 
The identity, as formulated in [20], is: 



If AS{f{ui,. . . ,Un-i)) = Eae5„_i(-l)'^/(Mi)'---'Mn-i))' ^"c/ (• • • )<o means 
keeping only non-positive powers of a Laurent polynomial in the variables U£, then 
the following equality holds: 



(4.24) < 



'n-l 



Yl (1 - ueum) AS MJ 2^+^ Yl (1 + "^""^ + ™™) 



l<e<m<n-l \l=l l<i<m<n-l 

/n— 1 \ n— 1 



<0 



AS n + ^ly-') = n n - ^^"')(^ + + 

\e=l / 1=1 l<l<m<n-l 



Applying to (g^S]), we find 



n-l 

Ui — Uj 



(4.25) Fn{t)= J] -^—^JY(l+tu^){l + TUi)\^-r^2.-. 

l<''<j<n-l «=1 

li t = 1, this coincides exactly with ()1.29p . What we have found is a very non-trivial combinatorial 
interpretation of this polynomial solution of gKZ, which goes beyond the special value q = e^'^*/^: 
for generic g, the sum of components of this solution reproduces the weighted enumeration of 
TSSCPPs with weight r = —q — q~^, proving a conjecture formulated in [22]. 

Next, we deduce from (|4.25p that the leading term of Fn{t) is 

n-l 

l<j<j<n-l * J i=l 

where the last equality is obtained by setting ui = and relabelling i ^ i — 1. With a bit more 
work, one can also prove that Fn{0) = t^~^ Fn-i{l / t) . 
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As a result, at r = 1, Fn{l), the sum of components, is An, and Fn{t)\in-i = Fn{0), the largest 
components, are A„_i. 

4.4.4. Refined enumeration. For general t, one can show that Fn{t) corresponds to a refined weighted 
enumeration of TSSCPPs, that is to adding extra boundary weights on the TSSCPPs. The details 
can be found in |28j . 

Remarkably enough, if one sets r = 1, then Fn{t) also coincides with the refined enumeration 
of ASMs. For ASMs, the refinement is simple to explain and consists in giving a weight of 
to an ASM whose 1 in the first row is at column i (which is essentially equivalent to keeping one 
spectral parameter free in the partition function of the six-vertex model with DWBC and sending 
the others to 1). This strange coincidence was observed in and proved via rather indirect 
arguments (making use of the Bethe Ansatz calculations of |91|). 

One can push this idea further: in 1986, Mills, Robbins and Rumsey conjectured that the double 
refinements of TSSCPPs and ASMs coincide [78]. Again, the double refinement of ASMs is easy to 
explain, consisting in weights t'^~^u^~^ for an ASM with a 1 at (1, i) and at (n, j) (which is equivalent 
to keeping two say horizontal spectral parameters free). The double refinement of TSSCPPs can be 
formulated in multiple ways and we refer the reader to appendix A of |34j for details. And in fact, 
using slight generalizations of the integrals above, one can give a direct proof of this conjecture and 
turn it into a theorem, as is done in |34j . 

5. Integrability and geometry 

The goal in this section is to reinterpret the families of polynomials that give rise to a solution 
of the gKZ equation in terms of geometric data. For simplicity we shall stick to cohomology (as 
opposed to -fC-theory), which means that we only consider rational solutions of the Yang-Baxter 
equation. In terms of the parameter g, it means taking an appropriate g — > ±1 limit (it is actually 
more convenient to consider q ^ —1 with our sign conventions). 

We shall consider three examples in what follows. The first one (matrix Schubert varieties) is 
there as warming up, though it already contains many of the key ingredients. The second example 
(orbital varieties) will correspond to the g — > — 1 limit of the solution of gKZ of section HI The third 
one (Brauer loop scheme) in fact contains the first two as special cases and is the most interesting 
one. But first we need to introduce some technology i.e. an algebraic analogue of equivariant 
cohomology for affine schemes with a linear group action: multidegrees [76]. Note that in this 
section we give almost no proofs and refer to the papers for details. 

5.1. Multidegrees. Our group will always be a torus T = (C^)^. T acts linearly on a com- 
plex vector space W . To a closed T-invariant sub-scheme X C we will assign a polynomial 
mdeg^y X G Sym(T*) = . . . , z^r] (here T* is viewed as a lattice inside the dual of the Lie 
algebra of T) called the multidegree of X. 

5.1.1. Definition by induction. This assignment can be computed inductively using the following 
properties (as in [38]): 

1. IfX = W = {0}, then mdeg^K X = 1. 

2. If the scheme X has top-dimensional components Xi, where nii > denotes the multiplicity 
of Xi in X, then mdeg^^ X = "^^mi mdegy^r X^. This lets one reduce from the case of 
schemes to the case of varieties (reduced irreducible schemes). 

3. Assume X is a variety, and H is a, T-invariant hyperplane in W. 

(a) UX ^ H, then mdegjy X = mdeg^(X n H). 
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(b) UX C H, then mdeg^ X = (mdeg^ X) ■ (the weight of T on W/H). 

One can readily see from these properties that mdeg^y X is homogeneous of degree codimvi/ X, and 
is a positive sum of products of the weights of T on W. 

5.1.2. Integral formula. We can also reformulate the multidegree as an integral (generalizing the 
idea that the degree of a projective variety is essentially its volume). To the torus T acting linearly 
on the complex vector space W is naturally associated a moment map /U, which is quadratic. With 
reasonable assumptions (i.e. that the multigrading associated to the torus action is positive), one 
can map generators of T* to complex numbers in such a way that this quadratic form is positive. 
Then one can formally write 

mdegi^ X = ^ ^ 

where it is implied that on both sides, this evaluation map has been applied. More explicitly, 
suppose that (2;j)j=i^,,,^„ is a set of coordinates on W which are eigenvectors of the torus action, 
with weights a^; then 

mdegvK X- ■'^ 



i=i 

where the must be evaluated to positive real numbers for the integral to make sense. So, the 
multidegree is just a Gaussian integral. The subtlety comes from the fact that in general X is 
a non-trivial variety to integrate on (in particular, it will be singular at the torus fixed point 0, 
the critical point of the function in the exponential!). In the case where X is a linear subspace of 
W , given by equations Xi^ = ■ ■ ■ = = 0, which is the only case where the integral is really an 
ordinary Gaussian integral, we compute immediately mdeg^y X = a^^ . . . ai^, as expected. 

5.2. Matrix Schubert varieties. Matrix Schubert varieties are a slight variation of the classical 
Schubert varieties, whose study started Schubert calculus. They are simpler objects to define, being 
afl&ne schemes (i.e. given by a set of polynomial equations in a vector space). 

5.2.1. Geometric description. Let g = g[(iV, C), b+ = {upper triangular matrices} C 0, b_ = 
{lower triangular matrices} C g, and similarly we can define the groups G = GL{N, C) and B± = 
{invertible upper/lower triangular matrices} C G. 

The matrix Schubert variety associated to a permutation a £ Sn is a closed sub-variety of g 
defined by the equations: 

(5.1) = {M e Q : rank Mj < ranker j, i,j = l,...,N} 

where M j is the sub-matrix above and to the left of {i,j), and o"^ is the transpose of the permutation 

i—i 

matrix of a. 

Explicitly, 5^- is defined by a set of polynomial equations of the form, determinant of a sub-matrix 
of M is zero, which express the rank conditions. 

One can also describe as a group orbit closure: 



(5.2) = B^a^B^ 



58 



Evidently, we have S„-i = Less evidently, the codimension of So- is the inversion number 

\a\ of a. 

Examples: 

• £7 = 1: then rankc^j = min(i, j) which is the maximum possible rank of a i x j matrix. So 

i I 

there is no rank condition and = Q- 

• 0" = (jQ, the longest permutation, (7o(i) = N-\-l — i: this time rank o""^ i = max{0,i + j — N) 



and we find that S^o is a linear subspace 

S^, = {M = (Mij-)i,i=i,...,7V : Mij = Vi,i, i + j<N} 



All the other matrix Schubert varieties are somewhere in between. For example, if a 
(4132), then 



a 



fO l\ 

10 

10 

\0 1 0/ 



s, 



(4132) 



M- 



11 



M- 



12 



M- 



13 



■ 1,...A ■ 







M21M32 - M22M31 = 



Note that among the multiple determinant vanishing conditions, one can extract in this 
case four that imply all the others. In general, there are simple graphical rules to determine 
which conditions to keep. 

There is a torus T = (C^)^^ acting on q, by multiplication on the left and on the right by 
diagonal matrices. The corresponding generators of its dual are denoted by xi,. . . ,X]y (left) and 
yi, . . . ,yN (right). In fact, since multiplication by a scalar does not see the distinction between 
left and right, the torus acting is really of dimension 2N — 1, and this amounts to saying that all 
multidegrees we shall consider only depend on differences Xj — yj. 

The are T-invariant; let us define the double Schubert polynomials Ho- to be their multidegrees: 



Note that E^~i{xi, ... ,XN\yi, ■■■ ,yN) 
play symmetric roles. 



= mdegg Sa 

r{yi, . . . ,yN\xi, . . . ,xn), so the two sets of variables 



5.2.2. Pipedreams. In [32j (see also [56])) a combinatorial formula for the calculation of double 
Schubert polynomials was provided. It can be described as configurations of a simple statistical 
lattice model coined in [56] (reduced) "pipedreams" . 

The pipedreams are made of plaquettes, similarly to loop models of section [3l The two allowed 



plaquettes are 



and 



(where the second one should be considered as two lines actually 



crossing). Furthermore, pipedreams have a specific geometry: they live in a right-angled triangle, 
as shown on the examples below. On the hypotenuse, the plaquettes are forced to be of the first 
type (and only one half of them is represented). 

A pipedream is reduced iff no two lines cross more than once. Alternatively, one can think of 
putting a non-local weight of zero to "bubbles" formed by two lines crossing twice. We assign 
weights to reduced pipedreams as follows: a crossing at row/column (i,j) has weight Xj — yj. We 
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can now state the formula for double Schubert polynomials: 
. . . ,a;Ar|yi, . . . ,yAr) = ^ 

reduced pipcdrcams of size N 
such that point i on vertical axis 
is sent to point cr{i) on horizontal axis 
i=l,...,N 



(weight of pipedream) 



A sketch of proof of this formula will be given in the next two sections. 
Examples: 

• £7 = 1: then only one pipedream contributes, the one with no crossings, of the type 

12 3 4 



■=1 



3 
4 

G = gq: again, only one pipedream contributes, this time with only crossings: 

12 3 4 















y 
























J 






















J 






J 







= n ~ 



Each factor of Xi — i/j has the meaning of weight of the equation Mij = 0. 
If o" = (4132), there are two (reduced) pipedreams: 




^(4132) = {xi - yi){xi - y2){xi - y3){x2 + X3 - yi - y2) 



Again, we recognize in each factor the weight of one equation of 5'(4i32) . In general, as long 
as Sa is a complete intersection (as many defining equations as the codimension), Hq- is just 
the product of weights of its equations (i.e. linear factors). 



5.2.3. The nil-Hecke algebra. In |33^ [32], double Schubert polynomials are related to a solution 
of the Yang-Baxter equation (YBE) based on the nil-Hecke algebra. We shall reformulate this 
connection here in terms of our exchange relation. 

Define the nil-Hecke algebra by generators ti, i = 1, . . . , N — 1, and relations 

j-2 



t 







^i^j — ijti \^ j\ ^ 1 
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As a vector space, the nil-Hecke algebra is isomorphic to C[5Ar] (the ti being identified with el- 
ementary transpositions). It has an obvious graphical depiction in which basis elements a £ Sn 
permute lines; the product is the usual product in Sn except that when lines cross twice, the result 
is zero: 

12 3 4 

12 3 4 



(2143)(1342) 




(2431) 



(1423)(3124) 



12 3 4 

(the expressions from right to left correspond to the pictures from bottom to top). The inversion 
number of a is graphically its number of crossings. 

Define 

Ri{u) = 1 + uti 

which is the associated solution to the rational Yang-Baxter equation, that is with additive spectral 
parameters: 

Ri{u)Ri+i{u + v)Ri{v) = Ri+i{v)Ri{u + v)Ri+i{u) 



Note that tj 




, and the identity is 1 




, that is 45 degrees rotated versions of our 



plaquettes, cf a similar remark in section 14.1.11 As usual, our solution of YBE also satisfies the 
unitarity equation: 

Ri{u)Ri{-u) = 1 

Now consider the formal object H = ^o-e5iv viewed as an element of nil-Hecke. We then 
claim that the following formulae hold: 

(5.3) H(xi, . . . ,XN\yi, ■ ■ ■ ,yN)Ri{xi+i - ^i) = E{xi, . . . ,Xi+i,Xi, . . . ,XN\yi, ■ ■ • ,y7v) 

(5.4) Ri{yi - yi+i)E{xi, . . .,XN\yi, ...,yN) = ^(xi, . . .,XN\yi, ■ ■ ■,yi+i,yi, ■ ■ ■,yN) 

which is nothing but exchange relations. Note that if we consider X as a function of the first set 
of variables xi, . . . , xj\f , the nil-Hecke algebra must act on the right on itself (ultimately, this is 
because of the transposition in (j5.ip and (|5.2p ). 

These formulae are obvious if one defines H in terms of pipedreams, that is if one writes S as [32] 



2/1 2/2 2/3 2/4 



xi 

X2 

X4 









J/ 






J/ 




J/ 


J/ 








61 



where each plaquette is an i?-matrix evaluated at yj — Xi (row i, column j), and each pipedream is 
interpreted as an element of nil-Hecke according to the convention of endpoints shown next to the 
equality. The proof of (j5.3p is then the usual argument (cf for example section I2.5.2p of applying 
the R matrix between rows i and i + 1, moving it across using YBE and checking that it disappears 
once it reaches the hypotenuse; and similarly for ()5.4p . 

To conform with the literature, we shall now focus on the first form (j5.3p . Writing it explicitly 
in components, we have to distinguish as usual cases. Let us denote as before Si the permutation 
of variables Xi and Xj+i. 

• If (t(z) < a{i + 1) (the lines starting at (z, i + do not cross): then a is not the image of tj 
acting by right multiplication, so that we find 

(5.5) = Si^cr 

i.e. Ho- is symmetric under exchange of Xi and Xj+i. 

• If a{i) > a{i + 1) (the lines starting at (i, i + 1) cross): this time we find that Ho- + (a^j+i — 
Xi)Ho' = SjHo, where a' is the unique preimage of a under right multiplication by ti. Thus, 

(5.6) Ho' = 9iHo 

where di = {si — 1), as before. Note that a' satisfies cr'(i) < cr'(f + 1), so that (j5.5p 

can be deduced from (|5.6p (the image of di is symmetric in Xj, Xj+i). 

In ()5.6p . a' has one fewer crossings than a. It is easy to show this way that, starting from Hoq = 
Y\i+j<N^^i ~ yj)^ where do is the permutation with the most crossings (the highest inversion 
number), one can compute any Ho by repeated use of ()5.6p . 

We now explain what the geometric meaning of (j5.6p is. This will prove that the pipedreams do 
compute the multidegrees of matrix Schubert varieties. 

5.2.4. The Bott-Samelson construction. Let Lj be the subgroup of G consisting of matrices which 
are the identity everywhere except in the entries Mat with a, 6 G {i, i + 1}, and B±^i = B± n Li. We 
use the following lemma, which is similar to lemma 1 of [58] (see also lemma 8 of [59]): (throughout 
the lemma the sign it is fixed) 

Let V be a left Li-module, and let X be a variety in V that is invariant under 
scaling and under the action of B±^i C Li. DeEne the map fi : (Lj x X)/B±^i V 
that sends classes of pairs {g, x) (where B±^i acts on the right on Li and on the left 
on X ) to g ■ X. If i^L is generically one-to-one, then 

mdegy Im fi = ^di mdegy X 

If on the other hand X is Li-invariant, then 

di mdegy X = 

Here we apply the lemma with V = g, Li acting on it by left multiplication, and X = S^- 
According to (|5.2p . X is i?„^j-invariant. Furthermore, if a{i) > a{i + 1), one can check that the 
map ^ is generically one-to-one, and the image is precisely S'o' where a = cr'ti. Hence ()5.6p holds. 
On the contrary, if a{i) < a{i + 1), X is simply Lj-invariant, and we conclude that ()5.5p holds. 

5.2.5. Factorial Schur functions. Pipedreams contain non-local information carried by the connec- 
tivity of the lines, just like the loop models of section [3l It is natural to wonder if there is a "vertex" 
representation similar to the six- vertex model for pipedreams. As this model is based on a Hecke 
algebra (which we have considered so far in the regular representation), there are in fact plenty 
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of vertex representations corresponding to various quotients of the Hecke algebra. The simplest 
one, analogous to the Temperley-Lieb quotient, is as follows. Suppose that cr is a Grassmannian 
permutation. By definition, this means that a has at most one descent, i.e. there is a A; such that 
i k implies a{i) < a{i + 1). Then one can group lines into two subsets depending on whether 
they start on the vertical axis at position i < k oi i > k. If we use different colors for them, we get 
a picture such as 

1 2 3 4 5 

1 



2 



3 



4 



5 

Several observations are in order. First, there is indeed no more non-local information in the sense 
that the connectivity of the endpoints can be entirely determined from the sequences of colors on 
the horizontal axis. Indeed, if the green endpoints are ai < • • • < ak and the red endpoints are 
6i < • • • < bn~k, there is only one Grassmannian permutation which is compatible with these colors, 
namely (ai, . . . , a^, 6i, . . . , hn-k)- Secondly, there is never any crossing below row fc, so does not 
depend on xt+i, ■ ■ ■ ,X]\f. According to ()5.5p . we also know that it is a symmetric polynomial of 
xi, . . . ,Xk- Thirdly, there are now five types of plaquettes: the four colorings of the non-crossing 
plaquette, but only one crossing plaquette (lines of the same color cannot cross each other!). In fact, 
by identifying red to occupied and green to empty, we recognize the six- vertex model configurations 
under the form of north-east going paths, in which the weight 62 = (red paths cannot straight 
to the right). Taking into accounts the weights, this is exactly the free fermionic five- vertex model 
considered at the end of section [2.4.31 (Fig. I13p . Furthermore, if we continue the red lines to the 
left so they all end on the same horizontal line, we find exactly the configurations that contribute 
to the Schur function sx, where A is encoded by the sequence of red and green endpoints at the top 
(which one can extend into an infinite sequence by filling with red dots at the left and green dots 
at the right). Finally, by comparing the weights, we find the equality: 

E„{xi, . . . , Xn\0, . . . , 0) = Sx{xi, ...,Xn) 

In this case, the ^^(xi, . . . , Xn\yi, ■ ■ ■ , Un) are usually called factorial Schur functions. They are 
essentially the same as double Schur functions, see [HH [80] for a detailed discussion. 

5.3. Orbital varieties. We shall move on to more sophisticated objects. In general, orbital vari- 
eties are the irreducible components of the intersection of a nilpotent orbit (by conjugation) with 
a Borel sub-algebra inside a Lie algebra. We cannot possibly reproduce the general theory of or- 
bital varieties, and here, we shall restrict ourselves to a very special type of orbital varieties which 
corresponds to the Temperley-Lieb algebra (for more general orbital varieties from an "integrable" 
point of view, see [M1E5])- 

5.3.1. Geometric description. We use the same notations as in section [5^21 but here all matrices 
are of even size = 2n. Furthermore, call n-|_ = {strict upper triangular matrices} C b+. 

Consider the affine scheme 

On = {M G n+ : = 0} 
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that is upper triangular matrices that square to zero. We have the fohowing description of its 
irreducible components: they are indexed by link patterns vr G P2n- To each vr we associate the 
upper triangular matrix 7r< with entries ('/r<)jj equal to 1 if i < j and paired by vr, otherwise. 
For example, 

/OOOOOOOIX 

/ ••. iooooo\ 


10 
10 


V •■■§/ 

Then the corresponding irreducible component Oj^ is given by the equations 

(5.7) = {M G g : = and rank Afl' < rankvr^i, j = 1, . . . , A^} 

where M j is the sub-matrix of M below and to the left of (i, j). Alternatively, O^^ can be defined 
as an orbit closure: 



TT = X — -\ \ 7r< 

2 3 4 5 6 



(5.8) Ott = • 7r< acts by conjugation 

Like all orbital varieties. On is equidimensional, and the codimension of in n+ is n{n — 1). 

There is a torus T = (C^)^"*"^ acting on n+, by conjugation by diagonal matrices and by scaling. 
The corresponding generators of its dual are denoted hy xi, . . . , xn and A. In fact, since conjugation 
by a scalar is trivial the torus acting is really of dimension N, and this amounts to saying that all 
multidegrees we shall consider only depend on differences Xi — Xj (and on a). 

Finally, the Oj^ are T-invariant; we define Qt^, the Joseph-Melnikov polynomial, to be 

il^ = mdeg„^ 

These are extended Joseph polynomials, in the sense that the original Joseph polynomials corre- 
spond to A = (no scaling action) . Melnikov is the one that initiated the study of these particular 
orbital varieties [75] . 

The specialization A = 1, = 0, corresponding to considering the action by scaling only, gives 
the degree of 0,r (i-e. the number of intersections with a generic linear subspace of complementary 
dimension; in the case of a complete intersection, it is just the product of degrees of defining 
equations). 

Examples: 

• If vr = 0", the base link pattern, then Oq" is a linear subspace where the upper-right n x n 
block is free while all other entries are zero. Thus, 



i.' 

I > i^o" — 

l<«<i<" n+l<i<j<2n 



(5.9) oo"=<|(q i)\ !^o"= n n (^+^* 



We note the similarity with (|4.10p . 
• A random example in size N = 6: 

(5.10) TT = /"^"ZI^^^ ^ 0^ = {m en+: ^ _ ^ 



Mi2 =M34 = M35 =M45 =0^ 
(M2)i6 = (M2)26 = 



1 2 3 4 5 6 

Again, there are simple rules to determine which equations are actually needed. Here we 
have six equations left which is the codimension, so that 

Ott = (a + Xi - X2){A + X3 - X4)(A + Xs - X5)(A + X4 - X5)(2A + Xi - Xq){2A + X2- Xq) 

Compare with the example at the end of section 14.2.11 
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5.3.2. The Temperley-Lieb algebra revisited. In general, Joseph polynomials are known to be related 
to representations of the corresponding Weyl group, in our case the symmetric group. Note that 
the latter is nothing but the limit g — > ±1 of the Hecke algebra. For the particular nilpotent orbit 
we have considered, we find unsurprisingly a representation of the Temperley-Lieb quotient of the 
symmetric group. However our extra scaling action makes the story much more interesting and 
produces a connection to the rational qKZ equation. 

Consider the Temperley-Lieb algebra with parameter r = 2 (or q = —1) acting as in section 
13.2.51 on the space of link patterns. The rational i?-matrix is of the form 

- {A-u) + uei 

Ri{u) = 





A + u 

where we recall that graphically, = and the identity is 1 = It can be deduced 

from the trigonometric i?-matrix (|4.ip by sending q — > — 1 in the following way 

The rational qKZ system is 

(5.11) Ri{xi - Xi+i)VLN{xi, . . .,xn) = 0.{xi, . . .,Xi+i,Xi, . . .,xn) 

(5.12) p~'^niy{xi, Xn) = {-lY~^VtN{x2, ■■■,XN,Xi+ 3a) 

It can be obtained from the (trigonometric) qKZ system (|4.2l4.3p by once again sending q to —1: 

The claim is that the vector of Joseph-Melnikov polynomials defined in the previous section, 
solves ()5.11l5.12p . and in fact coincides with the ^ ^ limit of (-1)"'^"'''^^ /'^ h~"'^"'~^^ ^ n , where 
is the solution of the trigonometric qKZ system that was discussed in section 14.21 

We can proceed analogously to section 14.2.11 and rewrite (|5.1ip in components by separating it 
into two cases: 

(5.13) {i,i + 1) not paired in vr : di ( ] = 



A -\- Xi Xi-^\ 

(5.14) (i, i + 1) paired in vr : — (a + Xj — Xi+i)di^Tr = 0^/ 



We shall now discuss the geometric meaning of (|5.13l5.14p . Since we know the base case ()5.9p . 
this will suffice to prove that Otv, the vector of multidegrees, satisfies the whole qKZ system. 
Still, it would be satisfactory to have a geometric interpretation of ()5.12p too. Unfortunately, it is 
currently unknown. Note that the effect of the r.h.s. of ()5.12p on multidegrees can be quite drastic: 
for example, starting from the example (|5.10p . one goes back to the base case 0^, cf ()5.9p : but so 
doing, two quadratic equations have turned into linear equations! 

Remark: one can write a rational qKZ equation which is analogous to (j4.5p but with additive 
spectral parameters. It should not be confused with the Knizhnik-Zamolodchikov (KZ) equation: 
the latter is recovered by the further limit A ^ 0, turning the difference equation into a differential 
equation. 
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5.3.3. The Hotta construction. The Hotta construction [37] is intended to explain the Joseph rep- 
resentation (Weyl group representation on Joseph polynomials) for an arbitrary orbital variety. 
When extended to our scaling action, it will produce the exchange relation (jS.lip . For more details 
see [251159]. 

The idea, as in the Schubert case, is to try to "sweep" with a subgroup Lj, here acting by 
conjugation, and to apply the lemma of section 15.2.41 We need to be a little careful that our 
embedding space, n^, is not Lj-invariant, so we must translate multidegrees in n+ to multidegrees 
in g. It is easy to see that 

mdegg = JJ(a + — Xj) mdeg^^ W W C n+ 

i>j 

so that the effect of sweeping with Li amounts to a "gauged" divided difference operator di'. 

1 / 

dif = — di{{A + Xi+i - Xi)f) = {A + Xi- Xi+i)di- 



A -\- X^-j-l Xi A -\- Xi 

Now, letting Li and -B+,i act by conjugation, start with an orbital variety (Dt^ which according 
to (|5.8p . is i3_|_^j-invariant. If one tries to sweep it with Li, one can in general produce non-upper 
triangular matrices. In fact, a small calculation shows that this occurs exactly if Mjj+i ^ 0. So, 
we must distinguish two cases: 

• If (i, i + 1) are not paired in vr, this means that the rank condition at {i, i -|- 1) in (j5.7p says 
that Mii^i = 0; which in turn implies that A-|-Xj — Xj+i | 0,^- Furthermore, in this case 

is Lj-invariant, so that = 0. This is exactly ()5.13p . 

• If {i, i + 1) is a pair in vr, then generically Ma^i 7^ in 0,^. We then proceed in two steps. 
Cutting: first we intersect O^r with the hyperplane Mjj+i = 0. Since the intersection is 
transverse, the effect on the multidegree is to multiply by the weight of the hyperplane, 
which is A + Xi — Xj+i. Sweeping: now we can sweep with Lj and we stay inside On. One 
can check that the map is generically one-to-one, and by dimension argument the image 
must be a union of orbital varieties (plus possibly some lower-dimensional pieces). The 
claim, which we shall not attempt to justify here, is that the orbital varieties thus obtained 
are exactly the proper preimages of vr under e^. So we find —di{{A + Xi — Xi+i)^},^) = 
Y,Tr'ytTT,ei-K'=TT^Tr'-, which is equivalent to ()5.14p . 



5.3.4. Recurrence relations and wheel condition. Several other constructions have simple geomet- 
ric meaning. We mention in passing here the meaning of recurrence relations and of the wheel 
condition. 

Set Xj+i = Xi + A. This gives the weight of to Ma^i. Roughly, this corresponds to looking at 
what happens when Ma^i — > 00 (this is more or less clear from the integral formula of section [5.1. 2 1 
for a more precise statement, see [57j). To remain inside On, writing the equations {M'^)ji^i = 
and {]VP)ij = 0, one concludes that one must have Mji = and Mj+ij = for all j. At this stage 
one notes that the entries Mji^i and Mij are unconstrained by the equations. Removing the i^^ 
and (i -|- 1)^* rows and columns reduces On to On-2- 

What we have just shown is that when Ma^i 00, being in On amounts to setting a certain 
number of entries to zeroes and then the result is some irrelevant linear space times On -2- One can 
be a bit more careful and do this reasoning at the level of each irreducible component: it is easy 
to see that only the components that have a pair (i,i + 1) survive (the others satisfy Ma^i = 0), 
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and that Oip^T^ gets sent to Ott- Finally, we get the recurrence relations: 
(5.15) . . = Xi + A, . . .) = < 



TT Im ip. 

i-l 2n 

TT(A + Xj-Xi) Y\ {2a + Xi-Xj) 
. . . ,Xi_i,a;i+2, . . . jXat) 



to be compared with ()4.14p . 

Similarly, if one sets x^ = Xj + 2a = Xi + A, i < j < k, then the equation (M^)jfc = cannot 
possibly be satisfied since it contains the infinite term MijMjk, so all multidegrees must vanish, 
which is nothing but the wheel condition. 

5.4. Brauer loop scheme. In this section we introduce a new affine scheme whose existence 
was suggested by the underlying integrable model. The subject has an interesting history, which 
we summarize now. In 2003, Knutson, in his study of the commuting variety, introduced the 
upper-upper scheme |55] ; one of its components is closely related to the commuting variety, and in 
particular has the same degree. In an a priori unrelated development, de Gier and Nienhuis [18] 
studied the Brauer loop model, a model of crossing loops which is completely similar to the one that 

we described in section [3.3.11 except crossing plaquettes ( ) are allowed. They found that the 

entries of the ground state (or steady state of the corresponding Markov process) are again integers, 
and observed empirically that certain entries coincide with the degrees of the irreducible components 
of the upper-upper scheme. This mysterious connection required an explanation. A partial one was 
given in [26], where the corresponding inhomogeneous model was introduced and it was suggested 
that this generalization corresponds to going over from degrees to multidegrees. Also, an idea of the 
geometric action of the Brauer algebra was given. But some entries remained unidentified. In [58], 
the Brauer loop scheme was introduced, and it was proved that its top-dimensional components 
produce all the entries of the ground state. Finally, in the recent preprint [59], it is shown that 
more generally, an appropriate polynomial solution of the gKZ system associated to the Brauer 
algebra produces multidegrees of these components. 

5.4.1. Geometric description. There are several equivalent descriptions of the Brauer loop scheme, 
see |58[ I59j. Here we try to present it in a way which best respects the underlying symmetries. 

Let N be an integer. Consider complex upper triangular matrices that are infinite in both 
directions and that are periodic by shift by (N,N): 

RzmodN = {M = {Mij)ij(zz : Mij = Mi+Nj+N G 

RzraoAN is an algebra. Let 5 = (5ij_i) G -Rz mod 7V be the shift operator. Then we define the 
algebra b+ to be the quotient 

b+ is finite-dimensional: dimb_(_ = A^^. Inside b+ we have its radical n+, which consists of (classes 
of) matrices with zero diagonal entries, and the group i?+ of its invertible elements, i.e. with 
non-zero diagonal entries. 

The Brauer loop scheme is then defined as 

En = {M en+: M"^ = 0} 
For a different point of view on the origin of E^, see section 1.3 of |59) . 
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In all that follows, we take N to be even, = 2n. We define a crossing link pattern to be an 
involution of Z/NZ without fixed points. Their set is denoted by P^n is of cardinality (2n — 1)!!. 
They have a similar graphical representation as link patterns, but in which crossings are allowed: 






per 

-'4 



The irreducible components of Ejy are known to be indexed by crossing link patterns |58^ I94j. If 
TT G P2n, then 

= {M e En : {M-^)i,i+N = {]S'P)j,j+N ^ i = j or i = 7r{j) mod N} 

Note that this definition (i) makes sense because (-/Vf^)j^j_|_7v, despite being undefined for a generic 
element of b+ (it is killed by the quotient by S^), is actually well-defined for elements of n+ and 
(ii) implies that the numbers (M^)j^j_|_7v always come in pairs for M G Ej\[, which is somewhat 
surprising (an elementary proof would be nice, as opposed to the proof of [5BJ ) • 

We can also describe E^^ as the closure of a union of orbits: 



7r{j) mod N} 



(5.16) E^ = B+-t7T< 

acts by conjugation, t7r< = {M G b+ : Mij 7^ : 
Finally, there is a conjectural description in terms of equations: 

(5.17) E^ = {M eb+:M^ = 0, {M\i+N = {M'^)^{i),-,{i)+N , rankM < rank7r<i-,i < i + N} 
where 7r< is any generic matrix in t7r<, for example the one made of zeroes and ones. 

Let us now turn to the torus action on b+. First there is the usual scaling action, with generator 
A. But the remaining A^-dimensional torus action is more difficult to explain. Let {xi)i^i be a set 
of formal variables satisfying the relations Xij^N = Xj + e, where e is another formal variable. Then 
the action of the full torus T = (C^)^"*"^ is defined by writing that 



wt(Mj. 



A + Xn 



i,j ez,i< j 



It is slightly non-trivial but true that E]\f and the E.,^ are T- invariant. We define the Brauer loop 
polynomials to be their multidegrees: 

= mdegf,^ E^ 

is a homogeneous polynomial in -|- 1 variables, which we can choose to be A, Xi — xi (i = 
2, . . . , A^) and e. Its degree is the codimension of E^^, which is 2n{n — 1). 

Examples: 

• The maximally crossed link pattern x"' corresponds to x"(0 = i + n. In this case, E^n is a 
linear subspace defined by the equations Mij = if j < i + n — 1. Thus, 



(5.18) 



X 




N i+n-1 



n n(-+ 

1=1 j=i+l 
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At the opposite end, we find the non-crossing link patterns of P2n- Among them, we have 
0" (or any of its rotations), which has ah pairings parahel. In the Brauer loop scheme, they 
play a special role: Knutson proved, in the context of the upper-upper scheme [55], that 
there is a Grobner degeneration of C„ x y to Eqk , where V is some irrelevant vector space, 
and Cn is the commuting variety: 

Cn = {{X, Y) G 0l(n, C):XY = YX} 

Explicit examples can become quite complicated; even in size = 4, we find, with the 
notation B = A — e: 



X 




^3 

E^2 = {M E n+ : M12 = M23 = M34 = M45 = 0} 

T^2 =(A + Xi- X2)(A + X2- X3)(A -I- X3 - X4)(B + Xi- Xi) 

M12 = M34 = 
M23M35 + M24M45 = 
M45M57 + M46M67 = 

M13M35-M24M46 = 0^ 
To2 =(a + Xi - X2)(a + X3 - X4) 

(A^ -|- AB -|- — BXi -|- AX2 + X1X2 — AX3 — X2X3 -|- BX4 — X1X4 -|- X3X4) 



0^ 




E^2 



M e n+ : 



M23 = M45 = 
M12M24 + M13M34 = 
M34M46 + M35M56 = 
M13M35-M24M46 = 0^ 



i<r — 

=^(0^) = f J = <J M G n+ : 

Tl2 =(A + X2 - X3)(B -I- X4 - Xi) 

(A^ -|- 2AB -|- BXi — AX2 — X1X2 + AX3 -|- X2X3 — BX4 -|- X1X4 — X3X4) 

The last two varieties are not complete intersections. 

5.4.2. The Brauer algebra. The Brauer algebra is defined by generators /j, ej, i = 1, 
relations 



A^-1, and 



3 



1 



(5.19) 

— djS-i fifj ~ fjfi ^ifj ~ fj^i K ~ j| ^ 1 

where indices take all values for which the identities make sense. Note that the fi are generators 
of the symmetric group Sn, whereas the generate a Temperley-Lieb algebra. 

There is an associated solution of the Yang-Baxter equation, namely 

a(a - u) + Auei + {1 - t/2)u{a - u)fi 



(5.20) 



Ri{u) 



(A + n)(A- (1-t/2)u) 



The Brauer algebra acts naturally on crossing link patterns. The graphical rules are the same 
as in the previous sections, and we shall not illustrate them again. If vr is viewed as an involution, 
then fi acts by conjugation by the transposition + 1), whereas e, creates new cycles + 1) 
and (7r(i),7r(z + 1)) - unless 7r(i) = i + 1, in which case it multiplies the state by r. 

We now claim that the vector T^v of multidegrees of the irreducible components of the Brauer 
loop scheme solves the qKZ system associated to the Brauer loop model [59], which we can write 

(5.21) Ri{xi - Xi+i)TN = SiTN yi £Z 

where Si permutes Xi+kN and Xi+i^kN for all k, on condition that the following identification of 
the parameters is made: 

(5.22) r = 

^ ^ 2a - e 

We could also add the cyclicity condition (which is obviously satisfied by T^v by rotational in- 
variance), but note that it could at most be marginally stronger than (j5.2ip (and in fact, it is 
not), because we have imposed (|5.21|) for all integer values of i. In other words, this is already an 
"affinized" version of the exchange relation because of our shifted periodicity property Xj+at = Xj+e. 

For future use, let us write in components ()5.2ip . We find the usual dichotomy: 

(5.23) TT{i)^i + l: -{A + Xi-X^+l){{A + B)di + S,)( ) = f.n 

\A + Xi-Xi+lJ 

(5.24) TT{i) = i + l: -{A + Xi-Xi+i){A + B + Xi+i-Xi)diT^ = {A + B) ^ T^/ 

with the convenient notation B = A — e. 

Remark: The Brauer algebra is the rational limit of the BWM algebra considered in |88j . 

5.4.3. Geometric action of the Brauer algebra. This is the most technical part of [59], which we 
shall only sketch here. Once again, it follows the same general idea of trying to "sweep" with a 
subgroup Li analogous to the Lj used so far. In our setting of infinite periodic matrices, this Lj 
consists of invertible matrices which are the identity everywhere except in the entries Mat with 
a,b £ {i + kN, i + 1 + kN} for some k; and B+^i = Lj H The additional subtlety comes from the 
fact that contrary to the case of orbital varieties, the condition = is not stable by conjugation 
by Li (i.e. (S"^) is not stable by conjugation by Lj, since the latter is outside b+). We shall have 
to reimpose one equation, namely {M'^)i^ii^N = 0, after sweeping. 

So, letting Li and B+^i act by conjugation, start with a component E,^ which according to (j5.16p . 
is S+^j-invariant. As usual we have to distinguish the two cases: 

• 7r(i) 7^ i + 1: in this case we can directly sweep with Li, and then cut with (M^)j_|_i^j_|_7v. As 
shown in [26l[58l[59] with varying levels of rigor and clarity, the result is precisely E-,^UEf..T^. 
So we find, at the level of multidegrees, 

-{A + B + Xi+i - Xi)diT^ = + T/.^ 

which reduces after a few manipulations to (j5.23p . 

• Tr{i) = 1 + 1: this time we first cut with Mj j+i = 0, sweep with Li (throwing away the 
Lj-invariant piece which cannot contribute to the multidegree calculation), and then cut 
again with (A/^)j_|_i j+Tv. We lost one dimension in the process, so the result cannot simply 
be a union of components. In fact, one can show that it is UttYt ei-7r'=7r ^ {iM'^)ii+N = 
(M2) 

i+ij+7V+i}- This results directly in the multidegree identity (j5.24p . 
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Recurrence relations can also be written and interpreted geometrically as in section 15.3.41 but 
we shall omit them for brevity. Let us simply mention the wheel condition: the T^^ satisfy 

^Tr{xi = X, Xj = X + A, Xk = X + 2 a) = i<j<k<i + N 

(cf the appendix B of ^8]). This stems geometrically from the equation {M'^)ik = (which is not 
killed by the quotient by since k < i + N), which cannot be satisfied when Mij,Mjk — > oo. 

5.4.4. The degenerate limit. There are two special values of the parameter r. One, on which we 
shall not dwell, is r = 1, which corresponds to e = in (j5.22p . This is the value for which T^r 
becomes the ground state eigenvector of an integrable transfer matrix, that of [18], and is the 
subject of [26l[58]. Various interesting results can be obtained, including sum rules. In particular, 
in the homogeneous limit — > multidegrees become degrees and we recover the integer numbers 
observed in [18]. 



However, there is another special point: r = 2. Indeed, we see that the coefficient of fi in 
the i?-matrix ()5.20p vanishes. In fact we recover this way the Temperley-Lieb solution of rational 
Yang-Baxter equation, which was used in connection with orbital varieties. What is the geometric 
meaning of this reduction? 

Plugging T = 2 into (j5.22p . we note that there is no solution for e, since e — > oo when r ^ 2. 
Let us carefully take the limit e — > oo in the polynomials T^r. We consider them as polynomials 
in zi, . . . , zn, a and e (this choice breaks the rotational invariance Zi — > Zj+i) and keep the former 
fixed while sending the latter to infinity. Geometrically, the situation is as follows. Let us compute 
the weights of the entries Mij of a matrix M S b+ in terms of zi, . . . , z^, A and B = A — e. Due to 
the periodicity, we can always assume i to be between 1 and A^, and due to the quotient we should 
then consider i < j < i + N . The result is that one finds two categories of entries: 



wt(Mj, 




- Zj l<i<j<N 

- Zj-N 1 < j - N < i < N 

This amounts to subdividing b+ (a vector space of dimension N"^) as 

f fx ^ \ 

b+ = 




V V — \J 

where U is an upper triangular matrix, and L a strict lower triangular matrix. These two pieces 
have quite different destinies as r ^ 2: the weights of U remain unchanged and are identical to 
those that we have used in section 15.31 for orbital varieties, whereas the weights of L diverge. It 
turns out that this limit corresponds to killing off this lower triangular part (as is clear from the 
definition of the multidegree as an integral, see section [5TL2]) ; and after taking the quotient, it is 
easy to see that we are left with b+, the algebra of upper triangular matrices. Let us call p this 
projection {U,L) i— > U from b+ to b+. 

At the level of components, here is what happens: if vr is a non-crossing link pattern, then 
p{Et^) = Ott) the corresponding orbital variety. Translated into multidegrees, this means: 

T^(zi, . . . , ZN,A, b) b"("~^) 0^(2:1, ...,zn,a) 
In general, crossing link patterns will project to lower dimensional i?+-orbit closures |59j . 

Interestingly, the matrix Schubert varieties and double Schubert polynomials of section 15.21 also 
appear as a special case in the r = 2 limit. Indeed, consider the permutation sector, that is the 
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subset of consisting of involutions such that vr({l, . . . , n}) = {n + 1, . . . , 2n}. Such involutions 
are in bijection with permutations a G 5^, according to = 7r(n + 1 — i) — n, i = 1, . . . , n: 




i — > vr 




flip 

i — > vr 




(the last flip is purely cosmetic and due to the fact that we have written labels of link patterns 
counterclockwise up to now). 

In the permutation sector, the structure of is as follows: 




If one discards the unconstrained entries (represented by one recognizes the pairs of n x n 
matrices {X,Y) in terms of which the upper-upper scheme of Knutson is defined [55]. In fact the 
union of E-,^ where vr runs over the permutation sector is exactly up to these irrelevant entries the 
upper-upper scheme. 

Now the projection of such components E.,^, that is here {X,Y) i-^ X, is essentially the corre- 
sponding matrix Schubert variety S^- There are various ways to see that; one can prove it rigorously 
by using the description in terms of orbits; or one can compare the defining equations ()5.ip of ma- 
trix Schubert varieties to the (conjectured) defining equations ()5.17p of Et^ (which, in the case of 
the permutation sector, reduce to the first conjecture of section 3 of [55]). In any case, the result 
is more precisely that there is a vertical flip, so that p(-E^) ^ SaCro- In terms of multidegrees, this 
means that 



T^(xi, . . . ,a;2n, a,b) ~ B 



B~*oo „n(n-l)-|(T| 



{A + Xi-Xj) JJ (A + Xi-Xj) 

l<i<j<" n+l<i<j<2n 

^cr (a -|- Xfn ■ ■ ■ 1 A -\- X\ \ Xfi-\-l ) • • • ) X2n) 



(note that the prefactor is nothing but r^o"). So the exchange relation (j5.3p satisfied by double 
Schubert polynomials (as well as the symmetric one (|5.4p ) should be somehow contained in the 
exchange relation (|5.2ip of Brauer loop polynomials. In order to see this, one must redefine the 
generator /j to tj = (1 — T/2)/j when taking the limit r — > 2. The operators ti then satisfy the 
nil-Hecke algebra, and one can check that in the permutation sector, where the part of the R- 
matrix never contributes when i ^ the exchange relation (|5.2ip reduces to ()5.3p if i < n or to 
([531) if i > n. 
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